To appear in The Astronomical Journal 

Preprint typeset using MTgX style emulateapj v. 5/2/1 1 



A SPECTROSCOPIC CATALOG OF THE BRIGHTEST (J < 9) M DWARFS IN THE NORTHERN SKY*'* 



m 
o 

X> 

PL, 
(N 



Oh: 
6 

(N 
> 

On 
On 

in 
o 

(N 



13 



SebastienLepine 1,2 , Eric J. Hilton 3 ' 4 , Andrew W. Mann 3 , Matthew Wilde 1 , Barbara Roias-Ayala 1 , KelleL. Cruz 1 - 2 ' 

& Eric Gaidos 4 

1 Department of Astrophysics, American Museum of Natural History, Central Park West at 79th Street, New York, NY 10024, USA; lepine@amnh.org, 

brojas-ayala@amnh.org, mwilde@amnh.org 
2 Department of Physics, Graduate Center, City University of New York, 365 Fifth Avenue New York, NY 10016, USA 
3 Institute for Astrophysics, University of Hawai'i, Honolulu, HI 96822, USA; hilton@ifa.hawaii.edu, amann@ifa.hawaii.edu 
4 Department of Geology & Geophysics, University of Hawai'i, 1680 East-West Road, Honolulu, HI 96822, USA; gaidos@hawaii.edu and 
5 Department of Physics and Astronomy, Hunter College, City University of New York, 695 Park Avenue, New York, NY 10065, kellecruz@gmail.com 

To appear in The Astronomical Journal 

ABSTRACT 

We present a spectroscopic catalog of the 1,564 brightest (J < 9) M dwarf candidates in the northern sky, 
as selected from the SUPERB LINK proper motion catalog. Observations confirm 1,408 of the candidates 
to be late-K and M dwarfs with spectral subtypes K7-M6. From the low (fj, >40 mas yr" 1 ) proper motion 
limit and high level of completeness of the SUPERBLINK catalog in that magnitude range, we estimate that 
our spectroscopic census most likely includes > 90% of all existing, northern-sky M dwarfs with apparent 
magnitude J <9. Only 682 stars in our sample are listed in the Third Catalog of Nearby Stars (CNS3); most 
others are relative unknowns and have spectroscopic data presented here for the first time. Spectral subtypes are 
assigned based on spectral index measurements of CaH and TiO molecular bands; a comparison of spectra from 
the same stars obtained at different observatories however reveals that spectral band index measurements are 
dependent on spectral resolution, spectrophotometric calibration, and other instrumental factors. As a result, 
we find that a consistent classification scheme requires that spectral indices be calibrated and corrected for each 
observatory /instrument used. After systematic correctio ns and a recalibration of the subtype-index relationships 
for the CaH2, CaH3, Ti05, and Ti06 spectral indices, we find that we can consistently and reliably classify 
all our stars to a half-subtype precision. The use of corrected spectral indices further requires us to recalibrate 
the £ parameter, a metallicity indicator based on the ratio of TiO and CaH optical bandheads. However, we 
find that our ( values are not sensistive enough to diagnose metallicity variations in dwarfs of subtypes M2 and 
earlier (±0.5dex accuracy) and are only marginally useful at later M3-M5 subtypes (±0.2dex accuracy). Fits of 
our spectra to the Phoenix atmospheric model grid are used to estimate effective temperatures. These suggest 
the existence of a plateau in the M1-M3 subtype range, in agreement with model fits of infrared spectra but at 
odds with photometric determinations of T e ff- Existing geometric parallax measurements are extracted from 
the literature for 624 stars, and are used to determine spectroscopic and photometric distances for all the other 
stars. Active dwarfs are identified from measurements of Ha equivalent widths, and we find a strong correlation 
between Ha emission in M dwarfs and detected X-ray emission from ROSAT and/or a large UV excess in the 
GALEX point source catalog. We combine proper motion data and photometric distances to evaluate the 
(U,V,W) distribution in velocity space, which is found to correlate tighly with the velocity distribution of G 
dwarfs in the Solar Neighborhood. However, active stars show a smaller dispersion in their space velocities, 
which is consistent with those stars being younger on average. Our catalog will be most useful to guide the 
selection of the best M dwarf targets for exoplanet searches, in particular those using high-precision radial 
velocity measurements. 

Subject headings: stars: low-mass, brown dwarfs - stars: late-type — surveys - catalogs - stars: fundamental 
parameters - stars: activity 



1. INTRODUCTION 

M dwarfs have become targets of choice for many exoplanet 
surveys. This is because low-mass planets (i.e. Earth- to 
Neptune-size) are easier to detect by the Doppler or transit 
techniques around stars of lower mass. The transits of smaller 
planets are also easier to detect when they occur in the also 
smaller M dwarfs. In addition, M dwarfs have much lower lu- 
minosities than the Sun, and their "habitable zones" (HZ) are 

* Based on observations collected at the MDM Observatory, located on 
Kitt Peak, and operated jointly by the University of Michigan, Dartmouth 
College, the Ohio State University, Columbia University, and the Univer- 
sity of Ohio. 

t Based on observations collected at the University of Hawaii 2.2- 
meter telescope, located on Mauna Kea, and operated by the University 
of Hawaii. 



closer in, which makes transits more likely to occur and ra- 
dial velocity variations easier to detect for planets in their HZ. 
Earth-like planets within the HZ of M dwarfs are thus emi- 
nently more detectable with current observationa l techniques 
than earth-like planets in the HZ of G dwarfs (Tart er et al.l 
120071: iGaidos et ai1l2007l) . M dwarfs are also the most plen- 
tiful class of stars, constituting the largest fraction (> 70%) 
of main se quence objects in the Galaxy and in the v icinity 
of the Sun (iReid. Gizis. & Hawlevl 120021 iCovev et alj|200a 
Boc hanskiet al.H2010l) . 

However, even nearby M dwarfs are generaly faint at the 
visible wavelengths where most planet searches are con- 
ducted, and most exoplanet detectio n techniques - w ith the 
notable exception of micro-lensing (Dong et al. 2009) - are 
currently restricted to relatively bright stars. This signifi- 
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cantly limits the number of M dwarfs that can be targeted 
in exoplanet surveys. Doppler searches in particular are usu- 
ally restricted to stars with visual magnitudes V < 12, and 
less than ~10% of late-K and early-M stars within 30 pc are 
currently being monitored by the lar ge-scale Doppler surveys 
dBufler et al.l200 6: Mayo r et al.l2009h . Ho wever, new survey s 
are pushing this limit to fainter magnitudes (lApps et al.l2010l) . 
and high-resolution spectrographs suitable for Doppler ob- 
servations at near-infrared wavelengths, where M dwarfs are 
relatively brighter, are being developed dTeradaet all 120081 : 
iBean et aUl2010t lOuirrenbach et al.ll2010t IWang et aUBOlOl) . 
In any case, only a fraction of all catalogued, nearby M dwarfs 
are bright enough to be included in radial-velocity monitoring 
programs. 

Tra nsit surveys, on the other hand, can include much fainter 



stars (Ir win et al 



120091) . However, because they have a 
much lower detection efficiency due to orbital inclination con- 
straints, they require extensive lists (thousands) of targets in 
order to detect any significant number of transit events. For 
transit surveys, the Solar Neighborhood census and its esti- 
mated «5,000 M dwarfs is therefore too small, and transit 
programs would greatly benefit from extending their target 
lists to much larger distance limits. 

A fundamental obstacle to progress has been the lack of 
a large, complete, and uniform catalog of bright M dwarfs 
suitable as targets for exoplanet programs. In particular, most 
catalogs and surveys of M dwarfs have focused on identifying 
the nearest objects, which are not necessarily the brightest. 
Whereas the Hipparcos catalog (Ivan Leeuwenll2007l) provides 
a near-complete census of solar-mass stars within 100 parsecs 
of the Sun, the bright magnitude limit of the catalog excludes 
all but the very nearest M dwarfs - although it lists stars as 
faint as V = 12-13, the Hipparcos catalog is complete only to 

y<8. 

The widely utilized Third C atalog of Nearby Stars, or 
CNS3 (lGliese& Jahreissl 119911) . which lists «3,800 stars, 
though predating the Hipparcos survey, has historically pro- 
vided a more complete list of M dwarf candidates in the So- 
lar Neighborhood. Many of the fainter stars in the CNS3 
have (grou nd-based) parallax measurements from a variety 
of sources (Ivan A ltena. L ee. & Hof fleit 1995). However, the 
CNS3 was largely compiled based on a photometri c anal- 
ysis th e high proper motion stars catalogued by iLuvten I 
dl979alfbT). and in larg e part using photometric d ata collected 
bv iGliesd (fl982l) and IWeisI <fl984l [1981 [19871) . The CNS3 
has been largely used in recent years to select M dwarf targets 
for exoplanet surveys dMarcv etal.l 12001 1 : iNaefetal.l 120031 
Butler et al.ll2004tlRiverall2005l: lEndl et al.ll2008b iBailev et alJ 
2009 ilAnglada- Escude et alj2012t lAnglada- Escude & Tuomil 



20121) . Unfortunately, the catalog suffers from various sources 
of incompleteness. These mainly consist of: (1) limited avail- 
ability of quality data for stars in the Luyten catalogs, at the 
time the CNS3 was compiled, (2) kinematic bias in the Luyten 
catalog due a relatively high (/i > 0.18" yr -1 ) proper motion 
limit at the low end, and (3) incompleteness of the Luyten 
catalogs even for stars with proper motions above the fiducial 
limit. 

Motivated mainly by the need to complete the census of 
the Solar Neighborhood, several surveys have since been 
conducted to identify the low-mass stars (mostly M dwarfs) 
suspected to be missing from the CNS3. These have in- 
cluded a re-analysis of the proper motion catalogs of lLuyteril 
( 1979b a) in light of high quality photometric data provided 



by the 2MASS survey dCutri et al.1 [20031) . This has led to 
the identification of hundreds of additional near by star can- 
didates that had prev iously been overlooked dReid & Cruzl 
120021: iReid et all 12004ft . In addition, cross-matching against 
2MASS and examination of Digitized Sky Survey images has 
uncovered significant (> 1') errors in many of the coordinates 
quoted in the Luyten catalo gs, which was hitherto preventin g 
efficient follow-up st udies iBakos. Sahu. & Nemethl (120021) : 
lSalim&Gouid1(l2003l) . 

Parallelling these efforts, 
veys have been conducted, 
proper motion stars missing 
with a focus on co mpleting 
Solar Neighborhood (|Lepine. Shara. & Ric 



new proper motion sur- 
mainly to find the high 
from the Luyten catalogs 
the stellar c ensus of the 
l200l 
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[ Deacon. Hamblv & Cooke|2005MLevinel l2005tlLe pine1l2005 
ISubasavage et al J 12005 allbb ILepind 120081) . In addition, some 
survey s have also been reaching to lower proper motion 
limits (ILepind I200H IReid. Cruz. & Alien! 120071: iBoyd eTall 
201 1), potentially extending the census of M dwarfs to larger 
distances. Recently, we have analyzed data from theSU- 
PERBLINK proper motion survey, which has a proper motion 
limit /i > 0.04" yr" 1 , with an emphasis on the identification 
of bright M dwarfs, rather than just nearby ones; our search 
has turned up 8,889 candidate M dwar fs with infrared mag- 
nitude J < 10 (Lepine & Gaidos 2011). Of these, we found 
that only 982 were previously listed in the Hipparcos catalog, 
and another 898 in the CNS3. Most of the other 7009 stars 
were not commonly known objects, and were identified as 
probable nearby M dwarfs for the first time. With its high 
estimated completeness , especially in the northern sky, the 
Lepine & Gaidql (1201 ll) census provides a solid basis for as- 
sembling an extensive and highly complete catalog of bright 
M dwarfs, suitable for exoplanet search programs. 

Not all M dwarfs, however, are equally suitable targets for 
planet searches. Some M dwarfs have significant photometric 
variability (flares, sp ots) which are affecting transit searches 
(lHartman et al.ll20~Tlb : some display chromospheric em ission 
affecting Doppler searches (llsaacson & Fischer] |20 10). Be- 
cause M dwarfs are relatively faint stars, they often require 
considerable investement of observing time on large tele- 
scopes to achieve exoplanet detection, and there is value 
in identifying subsets of M dwarfs that are intrinsically 
more likely to host detectable planets iHerrero et"aD (1201 ll) . 
In particular, one might be interested in selecting stars of 
higher metallicity which may harbor more massive planets 
dSousa et al J 120101) . or young stars with relatively luminous 
massive planets which would be easi er to detect through direct 
imaging (le.g. MugraueretaTll2010l) . In addition, one would 
like to avoid possible contaminants (e.g. background giants) 
or problematic systems (e.g. very active stars) in order to op- 
timize exoplanet survey efficiencies. 

Determining physical properties of the M dwarfs is also im- 
portant in order to better characterize the local populations of 
low-mass stars. This is especially true since proximity makes 
them brighter and thus more efficient targets for follow-up ob- 
servations and detailed study. Some of the bright M dwarfs 
may be close enough (d<20pc) to warrant inclusion in the 
parallax programs devoted to comple ting the census of low - 
mass stars in the Solar Neighborhood (le.g. Henry et al.l2 006). 
in which case it is also important that the candidates first be 
vetted through spectral typing. 

Spectral classification and analysis for a significant 
fraction of the low-mass stars in the CNS3 was per- 
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TABLE 1 

Survey stars: positions and photometry. 



Star name 


CNS3 a 


R.A.(ICRS) 


Decl.(ICRS) 




UDecI. 


Xray b 


hrl b 


FUV C 


NUV C 


V 


V 


J d 


H d 








(ICRS) 


(ICRS) 


" yr 1 


" yr 1 


cnts/s 




mag 


mag 


mag 


flag 


mag 


mag 


mag 


PM 100006+1829 




0.163528 


18.488850 


0.335 


0.195 








20.04 


11.28 


T 


8.44 


7.79 


7.64 


PM 1000 12+1 35 8S 




0.303578 


13.972055 


0.025 


0.144 








19.85 


11.12 


T 


8.36 


7.71 


7.53 


PM 100033+0441 




0.829182 


4.686940 


-0.024 


-0.085 








21.18 


12.04 


T 


8.83 


8.18 


7.98 


PM 100051+4547 


G12 


1.295195 


45.786587 


0.870 


-0.151 










9.95 


T 


6.70 


6.10 


5.85 


PM 100051+7406 




1.275512 


74.105217 


0.035 


-0.023 










10.63 


T 


7.75 


7.15 


6.97 


PM 100077+6022 




1.927582 


60.381760 


0.340 


-0.027 


0.1700 


-0.41 






14.26 


P 


8.91 


8.33 


8.05 


PM 100078+6736 




1.961682 


67.607124 


-0.045 


-0.091 










12.18 


P 


8.35 


7.72 


7.51 


PM 100081+4757 




2.026727 


47.950695 


-0.119 


0.003 


0.2190 


-0.27 


19.68 


18.91 


12.70 


P 


8.52 


8.00 


7.68 


PM 100084+1725 


GJ 3008 


2.113679 


17.424309 


-0.093 


-0.064 








19.24 


10.73 


T 


7.81 


7.16 


6.98 


PM 100088+2050 


GJ 3010 


2.224675 


20.840403 


-0.065 


-0.247 


0.0899 


-0.28 


21.07 


16.71 


13.90 


P 


8.87 


8.26 


8.01 


PM 100110+0512 




2.769255 


5.208822 


0.241 


0.061 






22.85 


20.58 


11.55 


T 


8.53 


7.88 


7.69 


PM 100113+5837 




2.841032 


58.617561 


0.056 


0.029 










11.21 


T 


8.02 


7.31 


7.13 


PM 100118+2259 




2.970996 


22.984573 


0.142 


-0.221 


0.4110 


0.28 




22.37 


13.09 


P 


8.86 


8.31 


7.99 


PM I00125+2142En 




3.139604 


21.713478 


0.189 


-0.290 










11.67 


T 


8.84 


8.28 


8.04 


PM 100131+7023 




3.298130 


70.398003 


0.045 


0.139 








19.93 


11.37 


T 


8.26 


7.59 


7.39 



a Designation in the Third Catalog of Nearby Stars I Gliese & Jahreiss 1 99 1 ) 

b X-ray flux and hardness ratio from the ROSAT all-sky points source catalog i Voges et al. 1999, 2000). 

c Far-UV and near-UV magnitudes in the GALEX fifth data release. 

d Infrared magnitudes from the Two Micron All-Sky Survey (Cutri et al.120031) . 



formed as part of the Palomar-MSU spectroscopic survey 
dReid. Hawlev. & Gizisl [T995L lHawlev. Gizis. & Reidl fl996h . 
hereafter PMSU. The survey has notably provided for- 
mal spectral classification for 1,971 of the fainter CNS3 
stars, confirming 1,648 of them to be nearby M dwarfs. 
More recent spectroscopic follow-up surveys have mainly 
focused on candidate nearby stars missing from the CNS3. 
Very high proper motion stars from th e Luyten catalogs 
dGizis & Reid 1 991 iReid & Gizisl 120051) . or stars discov- 
ered in the more recent proper motio n surveys | Scholz et all 
d2002b; iLepine. Rich. & Sharal d2003l) : IScholz et all d2005b : 
iRevle et alj d2006l ) have thus been targeted. Most notably, the 
"Meet the Cool Neighbours" program (hereafter MCN) has 
determined spectral subtypes for several hundred M dwarfs 
identified from the Luyten catalogs but not listed in the CNS3 
dCruz & Reidll2q02t ICruz"et al.ll2003t IReid et al.l l200l 120041: 
ICruz et al l 120071; IReid. Cruz. & Alien! 12007b . As with the 
other more recent surveys, the MCN program placed an em- 
phasis on the identification and classification of nearby, very- 
cool M dwarfs, most of which are however relatively faint 
and unsuitable for exoplanet surveys. It should be noted that 
while large numbers of M dwarfs have also been identified 
and classified as part of the spec troscopic follow-up program 
of the Sloan Dig ital Sky Survey (iBochanski et al.ll2005ll2010t 
IWest et aT]l201 ll) . most of them are relatively distant sources, 
and generally too faint for exoplanet surveys. 

This is why a signi ficant fraction of t he bri ght M dwarf can- 
didates published in Lepine & Gaidoi d201 lb had no available 
spectroscopic data at the time of release. In order to assem- 
ble a comprehensive database of M dwarfs targets suitable for 
exoplanet survey programs, we are now conducting a spectro- 
scopic follow-up survey of th e brightest M dwarf candidates 
from (Lep ine & Gaidosll201 lb . Our goal is to provide a uni- 
form catalog of spectroscopic measurements to confirm the M 
dwarf classification, and initiate detailed studies of their phys- 
ical properties, as well as the tayloring of exoplanet searches. 
In this paper, we present the first results of our survey, which 
provides data for the 1 ,564 brightest M dwarf candidates north 
of the celestial equator, with apparent near-infrared magni- 
tudes J<9. Observations are described in §2. Our spectral 
classification techniques are described in §3, and our model 



fitting and effective temperature determinations are given in 
§4. Metallicity measurements are presented in §5. Activity 
diagnostics are presented and analyzed in §6. A kinematic 
study informed by our metallicity and activity measurements 
is presented in §7, followed by discussion and conclusions in 
§8. 

2. SPECTROSCOPIC OBSERVATIONS 
2.1. Target selection 

Targets for the follow-up spectroscopic program were 
selected from the cata log of 8,889 bright M dwarfs of 
Lepine & Gaidoi (1201 lh . All stars are selected from the SU- 
PERBLINK catalog of stars with proper motions fi > 40 mas 
yr" 1 . Stars are identified as probable M dwarfs based on 
various color and reduced proper motions cuts; all selected 
candidates have, e.g., optical-to-infrared colors V — J > 2.7. 
The low proper motion limit of the SUPERBLINK catalog 
excludes nearly all background red giants. The low proper 
motion limit however also results in a kinematic bias, whose 
effects are discussed in §2.3 below. 

While some astrometric and photometric data have al- 
ready been compiled for all the stars, most lack formal spec- 
tral classification . Spec tral subtypes have been estimated in 
Lepine & Gaidoi (1201 lh only based on a calibrated relation- 
ship between M subtype and V — J color. However, the V 
magnitudes of many SUPERBLINK catalog stars are based 
on photographic measurements (from POSS-II plates); the re- 
sulting V—J colors have relatively low accuracy and are some- 
times unreliable. Besides from affecting spectral type esti- 
mates, unreliable colors can cause contamination of our sam- 
ple of M dwarf candidates by bluer G and K dwarfs, which 
would have otherwise failed the color cut. These are strong 
arguments for performing systematic spectroscopic follow-up 
observations, to provide reliable spectral typing and filter out 
G/K dwarfs (or any remaining M giant contaminants). 

A subsample of the brightest of the M dwarfs candidates, 
with apparent infrared magnitude J< 9 was assembled for 
the first phase of this survey. We also restricted the sam- 
ple to stars north of the celestial equator. This initial list 
contains a total of 1,564 candidates. All stars were in- 
discriminately targeted for follow-up observations, whether 
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FIG. 1 . — Examples of spectra collected at the MDM Obsevatory with either 
the Mklll or CCDS low-resolution spectrographs, and using the McGraw-Hill 
1.3-meter or Hiltner 2.4-meter telescope. All spectra covered the 6,200A- 
8,700A wavelength range with a resolution of 1.8A-2.4A per pixel, and a 
signal-to-noise ratio 20<S/N<30. 

or not they already had well-documented spectra. This 
would ensure completeness and uniformity, and allows com- 
parison of our sample with previous surveys. In particu- 
lar, o ur target list includes M dwarf c l assifi cation standards 
fromlKirkpatric k. Henry. & McCarthy (1991) which provide 
a solid reference for our spectral classification. The list in- 
cludes 557 stars that were previously observed in the PMSU 
spectroscopic survey, and 161 that were obseved and classi- 
fied as part of the MCN spectroscopic program (including 82 
stars observed in both the PMSU and MCN). 

Of the 1,564 M dwarf candidates, we found that 286 
had been observed at the MDM observatory by one of us 
(SL) prior to November 2008, as part of a separate spec- 
troscopic follow-up sur vey of very nearby (d<20pc) stars 
( Alpert & Lepinel |201 11) . The remaining targets were dis- 
tributed between our observing teams at the MDM Obser- 
vatory (hereafter MDM) and University of Hawaii 2.2-meter 
Telescope (hereafter UH), with the MDM team in charge of 
higher declination targets (<5 > 30) and the UH team in charge 
of the lower declination range (0 < S < 30). In the end we 
obtained spectra for all 1,564 stars from the initial target list. 

To check for any possible systematic differences arising 
from using different telescopes and instruments, we observed 
146 stars at both MDM and UH. We call this subset the "inter- 
observatory subset". Observations were obtained at different 
times at the two observatories. Data were processed in the 
same manner as the rest of the sample. 
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FIG. 2. — Examples of spectra collected at the University of Hawaii 2.2- 
meter telescope with the SNIFs spectrograph. Observations covered the 
5,200A-9,800A wavelength regime with a spectral resolution R~2000, and a 
signal-to-noise ratio 40<S/N<50. 

The full list of observed stars is presented in Table [TJ We 
used the standard SUPERBLINK catalog name as primary 
designation, however we also include the more widely used 
designations ( GJ, Gl, and Wo numbers ) for the 682 stars listed 
in the CNS3 dGliese & Jahreisslll991l) . The CNS3 stars are 
often well-studied objects, with abundant data from the liter- 
ature; a majority of them (580) were classified as part of the 
PSMU survey or MCN spectroscopic program. Another 56 
stars on our table which are not in the CNS3 were however 
classified as part of the PMSU survey or the MCN program. 
The remaining 821 stars are not in the CNS3, and were also 
not classified as part of the PMSU survey or MCN progra; 
these are new identifications for the most part, and little data 
existed about them until now. 

Table Q] lists coordinates and proper motion vectors for 
all the stars, along with astrometric parallaxes whenever 
available from the literature. The table also lists X-ray 
source counts from the ROSA T all-sky point source cat- 
alogs (IVoges et al.1 [1991 l2000b . far-UV (FUV) and near- 
UV (NUV) magnitudes from the GALEX fifth data re- 
lease, optical V magn itude from the SUPERBLINK catalog 
(Lepine & Shara 200 5b, and infr ared J, H, and K s magnitudes 
from 2MASS (ICutri etaDl2003h . More details on how thos e 
data were compiled can be found in Lepine & Gaidos (201 1). 
The optical (V band) magnitudes listed in Table Q]come from 
two sources with different levels of accuracy and reliability. 
For 919 stars in Table [TJ generally the brightest ones, the 
V magnitudes come from the Hipparcos and Tycho-2 cata- 
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TABLE 2 
Definition of Spectral Indices. 



150 



Index Numerator Denominator 



Reference 



CaH2 6814-6846 7042-7046 Reid, Hawlev. & Gizis (1995) 

CaH3 6960-6990 7042-7046 Reid. Hawl ev. & Gizis (1995) 

Ti05 7126-7135 7042-7046 Reid. Hawl ev. & Gizis (1995) 

VOl 7430-7470 7550-7570 Hawlev et al. (2002) 

Ti06 7550-7570 7745-7765 (Lepine, Rich, & Shara 2003) 

V02 7920-7960 8130-8150 (Lepine. Rich. & Shara 2003) 

logs. These are generally more reliable with typical errors 
smaller than ±0.1 magnitude; those stars are flagged "T" in 
Table Q] The other 645 objects have their V magnitudes esti- 
mated from POS S-I and/or POS S-II photographic magnitudes 
as prescribed in Lepine (2005). Photographic magnitudes of 
relatively bright stars often suffer from large errors at th e ~0.5 
magnitude level or more, in part due to photographic satura- 
tion; those stars are labeled "P" in Table Q] 

2.2. Observations 

Spectra were collected at the MDM observatory in a series 
of 22 observing runs scheduled between June, 2002 and April, 
2012. Most of the spectra were collected at the McGraw- 
Hill 1.3 -meter telescope, but a number were obtained at the 
neighboring Hiltner 2.4-meter telescope. Two different spec- 
trographs were used: the Mklll spectrograph, and the CCDS 
spectrograph. Both are facility instruments which provide 
low- to medium-resolution spectroscopy in the optical regime. 
Their operation at either 1.3 -meter or 2.4-meter telescopes is 
identical. Data were collected in slit spectroscopy mode, with 
an effective slit width of 1.0"to 1.5". The Mklll spectrograph 
was used with two different gratings: the 300 1/mm grating 
blazed at 8000A, providing a spectral resolution R~2000, 
and the 600 1/mm grating blazed at 5800A, which provides 
R~4000. The two gratings were used with either one of two 
thick-chip CCD cameras (Wilbur and Nellie) both having neg- 
ligible fringing in the red. Internal flats were used to cali- 
brate the CCD response. Arc lamp spectra of Ne, Ar, and Xe 
provided wavelength calibration, and were obtained for every 
pointing of the telescope to account for flexure in the sys- 
tem. The spectrop hotometric stand ard stars Feige 110, Feige 
66, and Feige 34 dOke et al.lll991t) were observed on a reg- 
ular basis to provide spectrophotometric calibration. Integra- 
tion times varied depending on seeing, telescope aperture, and 
target brightess, but were typically in the 30 seconds to 300 
seconds range. Between 25 and 55 stars were observed on a 
typical night. Spectra for a total of 901 bright M dwarf targets 
were collected at MDM. 

Additional spectra were obtaine d with the SuperN ova In- 
tegral Field Spectrograph (SNIFS; lLantz et al.ll2004l) on the 
University of Hawaii 2.2 m telescope on Mauna Kea between 
February 2009 and November 2012. SNIFS has separate but 
overlapping blue (3200-5600 A) and red (5200- 10000 A) spec- 
trograph channels, along with an imaging channel, mounted 
behind a common shutter. The spectral resolution is ^1000 
in the blue channel, and ^1300 in the red channel; the spatial 
resolution of the 225-lenslet array is 0.4". SNIFS operates 
in a semi-automated mode, acquiring acquisition images to 
center the target on the lenslet array, and bias images and cal- 
ibration lamp spectra before and after each target spectrum. 
Both twilight and dome flats were also obtained every night. 
Integration times depended on J magnitude but were 54 s for 
the faintest (J = 9) targets. Up to 75 target spectra were ob- 
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FIG. 3. — Histogram of the V — J color distribution of survey stars with 
spectral morphologies inconsistent with red dwarf of subtype K5 and later 
(in green). The distribution of V — J colors from the full sample is shown 
in red. Stars with spectra inconsistent with red dwarfs are close to the blue 
edge of the survey, which indicates that they are most probably contaminants 
of the color-magnitude selection, and not the result of mis-acquisition at the 
telescope. 

tained in one night. Spectra for 655 bright M dwarf targets 
were collected at UH with SNIFS. 

Spectroscopic data and results are summarized in Table [3] 
SUPEPvBLINK names are repeated in the first column and 
provide a means to match with the entries in Table Q] The 
second and third columns list the observatory and Julian date 
of the observations. The various spectroscopic measurements 
whose values are listed in the remaining columns are de- 
scribed in detail in the sections below. 

2.3. Reduction 

2.3.1. MDM data 

Spectral reduction of the MDM spectra was performed us- 
ing the CCDPROC and DOSLIT packages in IRAF. Reduc- 
tion included bias and flatfield correction, removal of the sky 
background, aperture extraction, and wavelength calibration. 
The spectra were also extinction-corrected and flux-calibrated 
based on the measurements obtained from the spectrophoto- 
metric standards. We did not attempt to remove telluric ab- 
sorption lines from the spectra. Many spectra were collected 
on humid nights or with light cirrus cover, which resulted in 
variable telluric features. However, telluric features generally 
do not affect standard spectral classification or the measure- 
ment of spectral band indices, since all the spectral indices 
and primary classification features avoid regions with telluric 
absorption. 

A more common problem at the MDM observatory was slit 
loss from atmospheric differential refraction. Although this 
problem could have been avoided by the use of a wider slit, 
the concomitant loss of spectral resolution was deemed more 
detrimental to our science goals. Instead, stars were observed 
as close to the meridian as observational constraints allowed. 
In some cases, however, stars were observed up to ±2 hours 
from the meridian, resulting in noticeable slit losses. Flue- 
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TABLE 3 

Survey stars: spectroscopic data 



Star name 


Observatory 


Julian Date 
2,450,000+ 


CaH2 r a 


CaH3 c 


Ti05 c 


voi c 


Ti06 r 


V02 r 


Sp.Ty. b 
index 


Sp.Ty. 

adopted 


c 


log g d 


Teff d 
K 


100006+1829 


MUM 


4791.75 




' 




' 








G/K 








PM 100012+1358S 


UH22 


5791.05 


0.706 


0.864 


0.788 


0.967 


0.944 


0.970 


0.14 


MO.O 


1.08 


5.0 


3790 


PM 100033+0441 


UH22 


5791.05 


0.580 


0.797 


0.679 


0.959 


0.888 


0.939 


1.38 


Ml. 5 


0.93 


4.5 


3510 


PM 100051+4547 


MDM 


5095.80 


0.603 


0.824 


0.664 


0.956 


0.911 


1.014 


1.10 


M1.0 


1.10 


4.5 


3560 


PM 100051+7406 


MDM 


5812.87 
















G/K 








PM 100077+6022 


MDM 


5838.82 


0.364 


0.620 


0.392 


0.905 


0.630 


0.798 


4.60 


M4.5 


0.90 


5.0 


3140 


PM 100078+6736 


MDM 


5099.83 


0.534 


0.789 


0.623 


0.905 


0.806 


0.929 


2.03 


M2.0 


0.96 


5.0 


3500 


PM 100081+4757 


MDM 


5098.84 


0.410 


0.700 


0.420 


0.871 


0.684 


0.814 


3.80 


M4.0 


1.01 


5.0 


3280 


PM 100084+1725 


UH22 


5791.06 


0.671 


0.842 


0.785 


0.970 


0.947 


0.970 


0.34 


M0.5 


0.91 


4.5 


3600 


PM 100088+2050 


MDM 


4413.72 


0.372 


0.646 


0.356 


0.893 


0.602 


0.793 


4.58 


M4.5 


0.99 


5.0 


3130 


PM 100110+0512 


UH22 


5050.03 


0.653 


0.839 


0.706 


0.962 


0.893 


0.925 


0.86 


M1.0 


1.16 


4.5 


3660 


PM 1001 13+5837 


MDM 


5811.90 
















G/K 








PM 100118+2259 


UH22 


5050.03 


0.439 


0.729 


0.424 


0.923 


0.694 


0.798 


3.50 


M3.5 


1.10 


4.5 


3260 


PMI0()125+2142En 


UH22 


5791.06 


0.721 


0.865 


0.851 


0.973 


0.967 


0.988 


-0.18 


MO.O 


0.80 


4.5 


3690 


PM 100131+7023 


MDM 


5812.89 


0.643 


0.816 


0.754 


0.973 


0.883 


1.038 


0.93 


M1.0 


0.88 


4.5 


3570 



a All spectral indices are corrected for instrumental effects, see Section 3.2. 

^ Numerical spectral subtype M evaluated from the corrected spectral band indices (not-i 

c Ha spectral index, comparable to equivalent width. 

^ Gravity and effective temperature from PHOENIX model fits. 

tuations in the seeing, which often exceeded the slit width, 
played a role as well. As a result, the spectrophotometric 
calibration was not always reliable, since the standards were 
only observed once per night to maximize survey efficiency. 
Flux recalibration was therefore performed using the follow- 
ing procedure. Spectral indices were measured and the spec- 
tra were classified using the classification method outlined be- 
low (§3.1). The spectra were then compared to classification 
templates assembled from M dwarf spectra from the Sloan 
Digital Sky Survey (SDSS) spectroscopic database. Each 
spectrum was divided by the classification template of the 
same alleged spectral subtype. In many cases, the quotients 
yielded a flat function, indicating that the spectrophotometric 
calibration was acceptable. Other quotients yielded residuals 
consistent with first or second order polynomials spanning the 
entire wavelength range, and indicating problems in the spec- 
trophotometric calibration. A second-order polynomial was 
fit through the quotient spectra, and the original spectrum was 
normalized by this fit, correcting for calibration errors. Spec- 
tral indices were then measured again, and the spectra reclas- 
sified; this yielded changes by 0.5 to 1.0 subtypes for w 20% 
of the stars (no changes for the rest). The re-normalization 
was then performed again using the revised spectral subtypes, 
and the procedure repeated until convergence for all stars. 

Finally, all the spectra were wavelength-shifted to the rest 
frames of their emitting stars. This was done by cross- 
correlating each spectrum with the SDSS template of the cor- 
responding spectral subtype. Spectral indices were again re- 
measured, and the stars re-classified. Any change in the spec- 
tral subtype prompted a repeat of the flux-recalibration pro- 
cedure outlined above, and the cross-correlation procedure 
was repeated using the revised spectral template until con- 
veregence. 

A sequence of MDM spectra is displayed in Figure[T] which 
illustrates the wavelength regime and typical quality. Note the 
telluric absorption features near 6,850A, 7,600A, and 8,200A. 
Signal-to-noise ratio is generally in the 30<S/N<50 range 
near 7500A. 



2.3.2. SNIFSdata 



SNIFS data processing was performed with the SNIF S data 
reduct ion p ipeline, which i s des cribed in detail in Baco n et ail 
(1200 11) and lAldering et all (120061) . After standard CCD pre- 
processing (dark, bias, and flat-field corrections), data were 
assembled into red and blue 3D data cubes. The data cubes 
were cleaned for cosmic rays and bad pixels, wavelength- 
calibrated using arc lamp exposures acquired immediately 
after the science exposures, and spectrospatially flat-fielded, 
using continuum lamp exposures obtained during the same 
night. The data cubes were then sky-subtracted, and the 
ID spectra were extracted using a semi-analytic PSF model. 
We applied corrections to the ID spectra for instrument re- 
sponse, airmass, and telluric lines based on observations of 
the Feige 66, Feige 1 10, BD+284211, or BD+174708 stan- 
dard stars (lOke et al.lll991l) . Because the SNIFS spectra are 
from an integral field spectrograph, operating without a slit, 
their spectrophotometry is significantly more reliable than the 
slit-spectra obtained at MDM. In fact, it is possible to perform 
synthetic photometry by convolving with the proper filter re- 
sponse. 

As with the MDM data, SNIFS spectra were shifted to 
the rest frames o f their emitting stars b y cross-correlation to 
SDSS templates (Bo chanski et al.l l2007) of the corresponding 
spectral subtype. Spectral indices were re-measured and the 
stars re-classified. This process was repeated if there was a 
change in the spectral subtype determination. 

A sequence of UH SNIFS spectra are displayed in Figure|2] 
which show the wavelength range and typical data quality. 
Signal-to-noise ratio is generally S/NwlOO near 7500A. 

3. SPECTRAL CLASSIFICATION 

3.1. Visual identification of contaminants 

We first examined all the spectra by eye to confirm the pres- 
ence of morphological features consistent with red dwarfs of 
spectral subtype K5 and later. The main diagnostic was the 
detection of broad CaH and TiO moleculars band near 7000A. 
Of the 1564 stars observed, 1408 were found to have clear 
evidence of CaH and TiO. The remaining 156 stars did not 
show those molecular features clearly, within the noise limit, 
and were therefore identified as probable contaminants in the 
target selection algorithm. 
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FIG. 4. — Comparison between our spectral band index measurements fo r a subset of 484 stars observ ed at MDM and UH, and the band index measurements 
for the same stars as reported in the Palomar-MSU spectroscopic survey of Reid, Hawlev. & Gizis 1 1995). Small but systematic offsets are observed, which are 
explained by differences in spectral resolution and spectrophotometric calibration between the different observatories. These offsets demonstrate that the spectral 
indices are instrument-dependent, but that the measurements can be corrected by observing large subsets of stars at the different obervatories. The red segments 
show the fits to the offsets, which are used to calibrate corrections for each observatory, here using the Palomar-MSU measurements as the standard. 



Most of these contaminants appeared to be early to mid- 
type K dwarfs, with a few looking like G dwarfs affected by 
interstellar reddening. A number of stars also displayed car- 
bon features consistent with low gravity objects, most prob- 
ably K giants. We suspect that many of the G and K dwarfs 
have inaccurate V-band magnitudes, making them appear red- 
der than they really are, which would explain their inclusion in 
our color-selected sample. Interstellar reddening would also 
explain the inclusion of more distant G dwarfs in our sample, 
due to their redder colors. An alternate explanation, however, 
might be that the targets were mis-acquired in the course of 
the survey, and that the spectra represent random field stars . 
Indeed the very large proper motion of the sources sometimes 
makes them difficult to identify at the telescope, as they of- 
ten have moved significantly from their positions on finder 
charts. Stars in crowded field are particularly susceptible to 
this effect. To verify this hypothesis, we compared the V — J 
color distribution of the contaminants to the distribution of the 
full survey sample (Figure [3] top panel); the fraction of con- 
taminant stars in each color bin is also shown (bottom panel). 
The two distributions are significantly different, with the con- 
taminants being dominated by relatively blue stars, and their 
fraction quickly drops as V—J increases. We can only con- 
clude that the contaminants are not mis-acquired stars, oth- 
erwise one would expect the two distributions to be statis- 
tically equivalent. Rather, the majority of the contaminants 
must have been properly acquired and are simply moderately 
red FGK stars that slipped into the sample in the photomet- 
ric/proper motion selection, as suggested first. 



We find an overall contamination rate of ~10% in our sur- 
vey, although most of the contamination occurs among stars 
with relatively blue colors. The contamination rate is ~26% 
for red dwarf candidates in the 2.7<V-7<3.0 color range, 
but this rate drops to ~8% for candidates with 3.0 < V—J < 
3.3, and becomes negligible (< 1%) in the redder 

V-J > 3.3 

candidates. The 156 stars identified as contaminants are in- 
cluded in Tableland Table[3]for completeness and future ver- 
ification. Spectroscopic measurements for these stars, such as 
band indices, subtypes, and effective temperatures, are how- 
ever left blank. 

3.2. Classification by spectral band indices 
3.2.1. M dwarf classification from molecular bandstrengths 

The spectra of M dwarfs are dominated by molecular bands 
from metal oxides (mainly TiO, VO), metal hydrides (CaH, 
CrH, FeH), and metal hydroxides (CaOH). The most promi- 
nent of these in the optical-red wavelength range (5000A- 
9,000A) are the bands from titanium oxide (TiO) and cal- 
cium hydride (CaH). The resulting opacities from those broad 
moleciular bands significantly affect the broad band colors 
and spectral energy distribution of M dwarfs dJonesj 119681 
lAllard et al.l200ttlKrawchuk, Dawson, & De RoberfislboOOft . 
Early atmospheric models of M dwarfs showed that the 
strength of the TiO and CaH bands depends on effective tem- 
perature, but also on surface gravity and metal abundances 
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FIG. 5. — Comparison between the corrected spectral band index measurements from MDM and the band index measurements of the same stars observed from 
UH. Observatory-specific corrections to the CaH2, CaH3, and Ti05 indices (see Figure |4} bring the MDM and UH values in close agreement. Small offsets 
are however observed for the Ti06, VOl, and V02 indices, which we could not calibrate against Palomar-MSU survey stars. The horizontal red lines show the 
adopted offsets which are used to correct the MDM values to bring them in line with the UH ones. Offsets are again believed to be due to inter-observatory 
differences in the spectral resolution and spectrophotometric calibration. 

ter quantify the classifi cation system, a number of "b and in- 
dices" were defined by iReid. Hawley. & Gizisl (119951) . which 
measure the ratio between on-band and off -band flux, for var- 
ious molecular bandheads. Calibration of these band indices 
against classification standards provide a means to objec- 
tively assign subtypes based on spectroscopic measurements. 
Originally, these band indices measured the strengths of var- 
ious features near 7,000, where the most prominent CaH 
and TiO features are found in early- type M dwarfs. These 
bands, however, become saturated in late-type M dwarfs, 
which makes their use problematic in later dwarfs. There 
is a VO band near 7,000A, located just between the main 
CaH and TiO features, which becomes prominent only at later 
subtypes; and band index measuring this feature was intro- 
duced as a primary diagnosis for the so-called "ultra-cool" M 
dwa rfs, and provided a classification scal e for subtypes M7- 
M9 (iKirkpatric k. Henry. & Simons! 1 19951) . Additional band 
indices associated with TiO and VO bands in the 7500A- 
9000A range, where molecular absorption develops at later 
subtypes, were introduced as s econdary classification features 
( Lepi ne. Rich. & Sha ra 2003). These form the basis of cur- 
rent classification methods based on optical-red spectroscopy. 

Nearby low-mass stars associated with the local halo pop- 
ulation have long been known to show peculiar banstrength 
ratios, and in particular to have weak T iO compared to CaH 
dMouldlll976l: iMould & McElrovlll978l) . A system to iden- 
tify and classify the metal-poor M subdwarfs based on the 
strength and ratio of CaH and TiO was introduced by iGizisi 



JMouldlll976l) . 

Molecular bands have historically been the defining diag- 
nostic and classification features of M dwarfs. For stars that 
have settled on the main sequence, one can assume that the 
surface gravity is entirely constrained by the mass and chem- 
ical composition. Leaving only the effective temperature and 
chemical abundances as general parameters in the classifica- 
tion and/or spectroscopic modeling. For local disk stars of 
solar metallicity, a classification system representing an ef- 
fective temperature sequence can thus be established based 
on molecular bandstrentghs. 

The detection and measurement of TiO and CaH molecular 
band s thus forms the basis for the M dwarf classification sys- 
tem dJoy & Abtlll974l) . Molecular bands become detectable 
starting at spectral subtype K5 and K7, the latest subtypes 
for K dwarfs (there are no K6, K8, or K9 subtypes). The 
increasing strength of the molecular bands then defines a se- 
quence running from M0 to M9. The strength of both TiO and 
CaH molecular bands reach a maximum around T e jj ~2700K. 
There is a turnaround in the correlation below this point, and 
molecular bands become p rogressively weaker at lower tem- 
peratures until they vanish (ICruz & Reidll2002h . The reversal 
and weakening is thought to be due to the condensation of 
molecules into dust, and their settling below the photosphere 
dJones&Tsuiilll997l) . 

Sequences of classification standard s were compiled in 
( Kirk patrick. Henry. & McCarthy| [T991). which identified the 
main molecular bands in the yellow-red spectral regime, 
where TiO and CaH bandheads are most prominent. To bet- 
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TABLE 4 

Coefficients of the spectral index corrections, by observatory. 



IDX 


A 1DX:PMSV 


BlDX.PMSU 


A/DXMDM 


BlDXMDM 


A]DX:UH 


BiDx.im 


CaH2 


1.00 


0.00 


0.95 


0.011 


0.92 


0.004 


CaH3 


1.00 


0.00 


0.90 


0.070 


1.00 


-0.028 


Ti05 


1.00 


0.00 


1.00 


0.000 


1.06 


-0.063 


VOl 






1.00 


0.040 


1.00 


0.000 


Ti06 






1.00 


-0.021 


1.00 


0.000 


V02 






1.00 


0.005 


1.00 


0.000 



d 19971). and expanded by ILepine. Shara. & Richl d2003l) and 
Lepin eTRich. & Sharal (120071) . as the spectroscopic census of 
M subdwarfs grew larger. In this system, the CaH band- 
strengths are used as a proxy of T e ff and determine spectral 
subtypes, while the TiO/CaH band ratio is used to evaluate 
metallicity. For that purpose, the ( Ti0 /CaH parameter, which 
is a function of the TiO/CaH band ratio, was introduced by 
Lepi ne. Rich. & Sharal (120071) as a possible proxy for metal- 
lic ity, and a tentative calibration with J Fe/Hl was presented 
by lWoolf. Lepine. & Wallerstein (2009). One of the main is- 
sues in the current M dwarf classification scheme, is that both 
TiO and CaH bandstrengths are used to determine the spectral 
subtype, whereas TiO is now believed to be quite sensitive to 
metallicity. This means, e.g., that moderately metal-rich M 
dwarfs may be assigned later subtypes than Solar-metallicity 
ones. In addition, the classification of young field M dwarfs 
may be affected by their lower surface gravities, which also 
tend to increase TiO bandstrengths and would thus yield to 
marginally later subtype assignments compared with older 
stars of the same T e jj. These caveats must be considered when 
one uses M dwarf spectral subtypes as a proxy for surface 
temperature. 

3.2.2. Definition and measurement of band indices 

The strength of the TiO, CaH, and VO molecular bands are 
measured using spectral band indices. These spectral indices 
measure the ratio between the flux in a section of the spec- 
trum affected by molecular opacity to the flux in a neighbor- 
ing section of the spectrum minimally affected by molecu- 
lar opacity. The latter section defines a pseudo-continuum of 
sorts, although M dwarf spectra do not have a continuum in 
the classical sense, because their spectral energy distribution 
strongly deviates from that of a blackbody, and is essentially 
shaped by atomic and molecular line opacities. 

We settle on a set of six spectral band indices: CaH2, CaH3, 
Ti05, Ti06, Ti07, VOl, an d VQ2. These band indixes 
which we previously used in ILepine. Rich. & Sharal (2003) 
to classify spectra collected at MDM, measure the strength 
of the most prominent bands of CaH, TiO, and VO in the 
6000 < A < 8500 regime. The spectral indices and are listed 
in Table [2] along with their definition. The CaH2, CaH3, 
and Ti05 indices are the same as tho se used in the Palomar- 
MSU survey, and were first defined in dReid. Hawley. & Gizisl 
19951). The TiQ6. VOl and V Q2 indices were introduced by 
Lepin e. Rich. & Sharal d2003l) to better classify late-type M 
dwarfs, whose CaH2 and Ti05 indices become saturated at 
cooler tempratuers and are not as effective for accurate spec- 
tral classification of late-type M dwarfs. Each spectral band 
index is calculated as the ratio of the flux in the spectral re- 
gion of interest (numerator) to the flux in the reference region 
(denominator), i.e.: 



Because the wavelength range for some indices is relatively 
narrow (especially the denominator for CaH2, CaH3, and 
Ti05) it is important that the spectra in which they are mea- 
sured have their wavelengths calibrated in the rest frame of 
the star, which is why special care was made to correct all 
spectra for any significant redshift/blueshift (see above). 

Because the measured molecular bandheads are relatively 
sharp, and because the spectral indices measuring them are 
defined over relatively narrow spectral ranges, the index val- 
ues are potentially dependent on the spectroscopic resolution, 
and may thus depend on the specific instrumental setup used 
for the observations. In addition, the index values may be 
affected by systematic errors in the spectrophotometric flux 
calibration, which can also be dependent on the instrument 
and/or observatory where the measurements were made. One 
way to verify these effects is to compare spectral index mea- 
surements of the same stars obtained at different observato- 
ries. Because of the significant overlap with the Palomar- 
MSU survey, we can use those stars as reference sample, and 
recalibrate the s pectral indices so that they are consistent to 
those reported in R eid. Hawlev. & Gizisl d 1995b . 

Our census have 557 stars in common with the PMSU spec- 
troscopic survey; these stars are all identified with a flag ("P") 
in the last column of Table [3] We identify 206 stars from the 
list observed at MDM and 28 1 stars from the list observed at 
UH which have spectral index measurements reported in the 
PMSU survey. The differences between our measured CaH2, 
CaH3, and Ti05 and those reported in the PSMU catalog are 
plotted in Figure [4] Trends and offsets confirm the existence 
of systematic errors, possibly due to differences in resolution 
and flux calibration. To verify the spectroscopic resolution 
hypothesis, we convolved the MDM spectra with a box kernel 
5-pixel wide; we found that indeed the MDM indices for the 
smoothed spectra had their offsets reduced by 0.01-0.02 units, 
bringing them more in line with the PMSU indices. We also 
observe that the MDM measurements tend to have a larger 
scatter than the UH ones; we suggest that this may be due to 
spectrophometric calibration issues with some of the MDM 
spectra, as discussed in §2.3.1. 

To achieve consistency in the measurements obtained at dif- 
ferent observatories, we adopt the values from the PMSU sur- 
vey as a standard of reference, and calculate corrections to the 
measurements from MDM and UH by fitting linear relation- 
ships to the residuals. A corrections to a spectral band index 
is thus applied following the general function: 



IDX C -Aiox-.obs IDXqbs + Bidx-.obs 



(2) 



f I(X)d\ 

ljy£ ~ num 



f, I(X)dX 



(i) 



where IDXobs represents the measured value of an index at 
the observatory OBS, and (Aidx-.obs^Bidx-.obs) are the coeffi- 
cients of the transformation from the observed value to the 
corrected one (IDX C ). Hence the corrected values of the in- 
dices CaH2, CaH3 and Ti05 for measurements done at MDM 
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FIG. 6. — Variation in the corrected spectral indices as a function of adopted spectral subtypes. The thin black lines show the adopted, revised calibrations for 
spectral classification. Blue circles represent the subset of classification standards from Kirkpatrick, Henrv. & McCarthy 1 1991) which we observed; the X-axis 
values of those data points are the formal spectral subtypes adopted in Kirkpatrick, Henry, & McCarthy ( 1991), while the Y-axis values are the spectral index 
measurements from our survey. For our own spectral typing scheme, we calculate the average of the subtype values for the CaH2 c , CaH3 c , Ti05 r , and Ti06 c 
indices given by the adopted relationships. The two VO indices have relatively weak leverage on early type stars due to the shallow slope in the relationships, 
and are not used for assigning spectral subtypes of our (mostly early-type) survey stars. 



are defined as: 

CaH2 c = AcaH2:MDMCaH2|yiDM + BcaH2:MDM (3) 
CaH3 c = AcaH3:MDMCaH3|yiDM + BcaH3:MDM (4) 
Ti05 c = ATi05:MDMTi05MDM + BTi05:MDM (5) 

The measurements from the PSMU survey are used as 
standards for these three indices, and we thus have by 

definition: Ac a H2:PMSU = AcaH3:PMSU = A T i05:PMSU = 1-0, 
#CaH2:PMSU = B Ca H3:PMSU = #Ti05:PMSU = 0.0. For OBS=MDM 

and OBS=UH, the adopted correction coefficients are listed in 
Table |4] The corresponding linear relationships are shown as 
red segments in Figure |4] 

To verify the consistency of the corrected spectral band in- 
dex values, we compare the corrected values for the stars ob- 
served at both MDM and UH (the inter-observatory subset). 
The differences are shown in Figure [5] We find the corrected 
values CaH2 c , CaH3 e , and Ti05 f to be in good agreement, 
with no significant offsets beyond what is expected from mea- 
surement errors. The corrected values of all three spectral in- 
dices are listed in Table [3] 

The Ti06, VOl, and V02 spectral index mesurements from 
MDM and UH are also compared in Figure [5] Those were 
not measured in the PMSU survey, since the indices were 
introduced later (Lepine, Rich, & Shara 2003), and thus are 
displayed here without any correction. Small but significant 
offsets between the MDM and UH values again suggest sys- 
tematic errors due to differences in spectral resolution and flux 



calibration. This time we adopt the UH measurements as fidu- 
cials, and determine corrections to be applied to the MDM 
data. The corrections are listed in Table |4]and the correspond- 
ing linear relationships are displayed in Figure |5] as red seg- 
ments. The corrected values of the three spectral indices are 
also listed in Table [3] 

The scatter between the MDM and UH values, after correc- 
tion, as well as the scatter between the UH/MDM and PMSU 
values, provide an estimate of the measurment accuracy for 
these spectral indices. Excluding a few outliers, the mean 
scatter is «0.02 units (ler) for the CaH2, CaH3, Ti05, Ti06, 
and VOl indices, and ?»0. 04 units for V02. This assumes that 
M stars do not show any significant changes in their spectral 
morphology over time, and that the spectral indices should 
thus not be variable. 

3.2.3. Spectral subtype assignments for K/M dwarfs 

Because of the correlation between spectral subtype and the 
depth of the molecular bands, it is possible to use the values of 
the spectral band indices to estimate spectral subtypes. This 
only requires a calibration of the relationship between spec- 
tral index values and the spectral subtypes, in a set of stars 
which were classified by other means, e.g., classification stan- 
dards. The system adopted in this paper uses the spectral in- 
dic es listed in T able [2] a nd follows the methodolo gy outlined 
in (iGizisll 19971) and (Lep ine. Rich. & Sharall2003l) . Relation- 
ships are calibrated for each spectral index, and spectral sub- 
types are calculated from the mean values obtained from all 
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relevant/available spectral indices. The mean values are then 
be rounded to the nearest half integer, to provide formal sub- 
typing with half-integer resolution. The system is extended to 
late-K dwarfs as well: an "M subtype" with a value < 0.0 sig- 
nifies that star is a late-K dwarf: the star is classified as K7 for 
an index value w — 1.0 and as K5 for an index value « -2.0 
(note: there is no K6 subtype for dwarf stars, and K7 is the 
subtype immediately preceding MO). 

The original spectral-index classification method for M 
dwarfs/subdwarfs is based on a relationship between sub- 
type and with the CaH2 index, which measures one of the 
most prominent band at all spectral subtypes, and notabl y dis- 
plays the deepest bandhead in meta l-poor M subdwarfs lGizisI 
( 119971) : ILepine. Rich. & Sharal d2003h . The original relation- 
ship is: [SpTy] CaH2 = 10.71 -20.63 CaH2 + 7.91 (CaH2) 2 . To 
verify this relationship, we estimated spectral subtypes from 
our corrected indice s CaH2 r for 16 spectroscopic calibra - 
tion standards from Kirkpatri ck. Henry. & McCarthyl (119911) . 
which were observed as part of our survey, and span a range of 
spectral subtypes from K7.0 to M6.0. We found small but sig- 
nificant differences in our es timated spectral subtypes and the 
values formally assigned by Kirkpatrick, Henry, & McCarthy 
(119911) : subtypes estimated from the iGizisI (119971) relationship 
tend to systematically underestimate the standard subtypes by 
ps 0.5 units for stars later than M3. To improve on the index 
classification method, we performed a x 2 polynomial fit to 
recalibrate the relationship, obtaining: 

[SpTy] CaH2 = 11.50-21.71 CaH2 c + 7. 99 (CaH2) 2 (6) 

which does correct for the observed offsets at later types. Us- 
ing this this relatioship as a starting point, and guided by the 
formal spectral subtype from the classification standards, we 
performed additional \ 2 polynomial fits to calibrate an index- 
subtype relationships for CaH3: 

[SpTy] CaH3 = 18.80-21.68 CaH3 c (7) 

Where the corrected values of the spectral bands indices (see 
Eqs.2-5) are used. The relation ship s are slightly different 
from t hose quoted in IGizisI (119971) and ILepine. Rich. & Sharal 
(2003) but are internally consistent to each other, whereas an 
application of the older relationships to our corrected band in- 
dex measurements would yield internal inconsistencies, with 
subtype difference up to 1 spectral subtype between the rela- 
tionships. 

The ratio of oxides (TiO, VO) to hydrides (CaH,CrH,FeH) 
in M dwarfs is known to vary significan tly with metallicity 
(lGizislll997t ILepine. Rich. & Sharall2007l) . In the metal-poor 
M subdwarfs, it is the oxides bands that appear to be weaker, 
while hydride bands remain relatively strong (in the most 
metal-poor ultrasubdwarfs, or usdM, the TiO bands are al- 
most undetectable). Therefore it makes sense to rely more on 
the CaH band as the primary subtype/temperature calibrator. 
The same [SpTy]c a H2 and [SpTy]c a H3 relationships should be 
used to determine spectral subtypes at all metallicity classes 
(i.e. in M subdwarfs as well as in M dwarfs). 

Because the TiO and VO bands are also strong in the metal- 
rich M dwarfs, it is still useful to include these bands as sec- 
ondary indicators, to refine the spectral classification. In the 
late-type M dwarfs, in fact, the CaH bandheads are saturating, 
and one has to rely on the TiO and VO bands. In fact, the VO 
bands were originally used to diagnose and calibrate ultracool 
M dw arfs of subtypes M7-M9 (Kirkpatrick, Henry, & Simons 
1995). The main caveat in using the oxide bands for spectral 
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classification is that this can potentially introduce a metallic- 
ity dependence on the estimated spectral subtype, with more 
metal-rich stars being assigned later subtypes than what they 
would have based on the strength of their CaH bands alone. 
In any case, because our sample appears to be dominated 
by near-solar metallicity stars, we calibrate additional rela- 
tionships between subtype and the Ti05 and Ti06 bands in- 
dices. We first recalculate the subtypes by averaging the val- 
ues of [SpTy] CaH2 and [SpTy] CaH2 , and perform a \ 2 fit of 
the Ti05 and Ti06 indices to the mean subtypews calculated 
from CaH2 and CaH3, finding: 

[SpTy] T . 05 = 7.83 - 9.55 Ti05 c (8) 
[SpTy] T . 06 = 9.92-15.68 Ti06 c + 21.23 (Ti06 c ) 2 

- 16.65 (Ti06 c ) 3 (9) 

where again the corrected band indices are used. The rel- 
atively sharp non-linear deviation in the Ti06 distribution 
around M3 forces the use of a third order polynomial in the 
fit. 

We also determine the relationships for the VOl and V02 
band indices. This this after recalculating the subtypes from 
the average of [SpTy] CaH2 and [SpTy] CaH2 , [SpTy] 

Ti05' aI1C ^ 

[SpTy] T . 06 , a x 2 fit again yields: 

[SpTy] VQ1 = 69.8-71.4 [V01] c (10) 
[SpTy] VQ2 =9.56-12.47 [V02] c + 22.33 ([V02] c ) 2 

-19.59 «V02] C ) 3 . (11) 

The VO indices however make relatively poor estimators of 
spectral subtypes for our sample, mainly because the shallow 
slope at earlier subtypes provides little leverage. The V02 in- 
dex also shows unexpectedly large scatter in the MDM spec- 
tra, including in the classification standard stars, which we 
suspect is due the fact that the index is defined very close to 
the red edge of the MDM spectral range and is thus more sub- 
ject to statistical noise and flux calibration errors. We there- 
fore do not include [SpTy] vol2 and [SpTy] V02 in the final de- 
termination of the spectral subtypes. 

Fig.|6]plots all the corrected spectral band indices as a func- 
tion of the adopted spectral subtype. The relatively small scat- 
ter (ra 0.02) in the distribution of [CaH2] c , [CaH3] c , [Ti05] c 
and [Ti06] c demonstrate the internal consistency of the spec- 
tral type calibration for the four indices. All four relationships 
have an average slope « 10, which means that since those 
indices have an estimated measurement accuracy of w 0.02 
units, the spectral subtypes calculated by combining the four 
indices should be accurate to about ±0.10 subtype assuming 
that the measurement errors in the four indices are uncorre- 
cted. While this would make it possible to classify the stars 
to within a tenth of a subtype, we prefer to follow the general 
convention and assign spectral subtypes to the nearest half in- 
teger. 

To verify the consistency of the spectral classification, we 
compare the spectral types evaluated independently for the 
list of 141 stars observed at both MDM and UH. We find 
that 82% of the stars end up with the same spectral type as- 
sigments from both observatories, i.e. they have spectra as- 
signed to the same half-subtype. All the other stars have 
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TABLE 5 

Distribution by spectral subtype of the 
1564 survey stars 



Spectral subtype 


N 


G/K a 


160 


K7.0 


27 


K7.5 


101 


MO.O 


177 


M0.5 


160 


M1.0 


152 


M1.5 


147 


M2.0 


141 


M2.5 


119 


M3.0 


125 


M3.5 


125 


M4.0 


72 


M4.5 


36 


M5.0 


13 


M5.5 


4 


M6.0 


2 


M6.5 


2 


M7.0 


1 



Stars identified as earlier than K7.0 and/or with no de- 
tected molecular bands. 

classifications within 0.5 subtypes. This is statistically con- 
sistent with a lcr error of ±0.18 on the spectral type deter- 
mination, slightly larger than the assumed 0.1 subtype preci- 
sion estimated above. This suggests a 3er error of about ±0.5, 
which justifies the more conservative use of half-subtypes as 
the smallest unit for our classification. 

The resulting classifications based on the CaH and TiO 
band index measurements are listed in Table [3] The numer- 
ical spectral subtype measured from the average of the band 
indices is listed to 2 decimal figures. These values are rounded 
to the nearest half integers to provide our more formal spectral 
classifications to a half-subtype precision. The non-rounded 
values are however useful for comparison with other physical 
parameters as they provide a continuous range of fractional 
values; these fractional subtype values are used in the analy- 
sis throughout the paper 

A histogram of the distribution of spectral subtypes is 
shown in Figure [7] with the final tally compiled in Table [5] 
Most of the stars in our survey have subtypes in the M0.0- 
M3.0 range. The sharp drop for stars of subtypes K7.5 
and K7.0 is exp l ained by the color selection used in the 
Lepine & Gaidos (201 1) catalog (V-J > 2.7) which was orig- 
inally intended to select only M dwarfs; note that stars with 
subtypes earlier than K7.0 are also excluded from the graph, 
and are probably contaminants of the color selection in any 
case. Our deficit of K7.0 and K7.5 stars however spectra 
demonstrates that the adopted selection criterion is efficient 
in excluding K dwarfs from the catalog. The distribution 
of spectral subtypes also shows a marked drop in numbers 
for subtypes M4 and later. This is a consequence of the 
relatively bright magnitude limit (J < 9) of our subsample, 
combined with the low absolute magnitudes of late-type M 
dwarfs, which excludes most late-type stars from our survey, 
since these tend to be fainter than our magnitude limit, even 
relatively nearby ones. 

3.3. Semi-automated classification using THE HAMMER 

To verify the accuracy and consistency of spectral typing 
based the spectral-index method described above, we per- 
formed in dependent spectra l classification using the Ham- 
mer code dCovev et al .112007b . The Hammer was designed to 
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subtype M 

FIG . 7 . — Distribution of spectral subtype M, for the stars in our survey, with 
K5=-2 and K7=-l. Most stars are found to be early-M objects. The sharp drop 
at subtypes K7.5 and earlier is a consequence of our initial V -J > 2.7 color 
cut, which was meant to select only stars cool enough and red enough to be 
M dwarf. The sharp drop for subtypes M4 and later is a consequence of the 
high magnitude limit (J < 9) of our survey, which restricts the distance range 
over which late-M stars are selected. The magnitude limit also explains the 
slow drop in numbers from subtypes M0 to M3, whereas one would normally 
expect the lower-mass M3 stars to be more common than earlier-type objects 
in a volume limited sample. 
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FIG. 8. — Comparison of spectral types assigned by the spectral index 
method and those assigned by the Hammer code, which is used for stellar 
classification in the Sloan Digital Sky Survey. Subtypes generally agree to 
within the advertized precision of the Hammer (±1 subtype, illustrated by 
the slanted lines), although the Hammer subtypes are marginally later than 
the spectral index subtypes, by 0.27 subtype on average. 

classify stars in the Sloan Digi tal Sky Survey Sp ectroscopic 
database, including M dwarfs dWest et al.ll201 ll) . The code 
works by calculating a variety of spectral-type sensitive band 
indices, and uses a best-fit algorithm to identify the spec- 
tral subtype providing the best match to those band indices. 
For late-K and M dwarfs, spectral subtypes are determined to 
within integer value (K5, K7, M0, Ml, M9). 

However, to ensure that the automatically determined spec- 
tral types were accurate we used the manual "eye check" 
mode of the Hammer (version 1.2.5). This mode is typically 
used to verify that there are no incorrectly typed interlopers. 
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The Hammer allows the user to compare spectra to a suite 
of tem plate spectra to determine the best match. IWest et all 
d201 ll) have found that for late-type M dwarfs, the automatic 
classifications were systematically one subtype earlier than 
those determined visually. Our analysis confirms this offset, 
and we therefore disregarded the automatically determined 
Hammer values to adopt the visually determined subtypes. 

The resulting subtypes are listed in Table [3] Some 170 stars 
were not found to be good fits to any of the K5, K7, or M type 
templates, and thus identified as early-K or G dwarfs. This 
subset includes all of the 156 stars that were visually identified 
as non-M dwarfs on first inspection (see §3.1). The remaining 
14 stars were initially found to be consistent with late-K stars, 
and classified as K7.0 and K7.5 objects using the spectral in- 
dex method described in §3.2 above; we investigated further 
to determine why the stars were classified as early K using 
the Hammer. On closer inspection, we found that 3 of the 
stars are indeed more consistent with mid-K dwarfs that K7.0 
or K7.5, and we thus overran the spectral index classification 
and reclassified them as "G/K" in Table [3] For the other 1 1 
stars flagged as mid-K type with the Hammer, we determined 
that the stars do show significant evidence for TiO absorption, 
which warrants that the stars retain their spectral-index clas- 
sifiction ofK7.0/K7.5. 

In the end Table[3]lists 159 stars from our initial sample that 
are identified and early-type G and K dwarf contaminants, and 
most likely made our target list due to inaccurate or unreliable 
V — J colors. The remaining 1405 stars are formally classified 
as late-K and M dwarfs. 

A comparison of spectral subtypes determined from the 
spectral-index and Hammer methods is shown in Figure [8] 
Because the Hammer yields only integer subtypes, we have 
added random values in the -0.4,0.4 range to facilitate the 
comparion. Slanted lines in Figure |8]show the range expected 
if the two classification methods (spectral index, Hammer) 
agree to within 1.0 subtypes. There is however a mean off- 
set of 0.26 subtypes between the spectral index and Hammer 
classifications, with the Hammer subtypes being on average 
slighlty later. 

Figure [8] also reveals a number of outliers with large dif- 
ferences in spectral subtypes between the two methods. We 
found 48 stars with differences in spectral subtyping larger 
than ±1.5. The spectra from these stars were examined by 
eye: except for one star, we found the band-index classifica- 
tions to agree much better with the observed spectra than the 
Hammer-determined subtypes. The one exception is the star 
PM 111055+4331 (Gl 412B) which the band index measure- 
ments classify as M6.5; in this case the Hammer determined 
subtype of M 5.0 appears to be more accurate. This excep- 
tion likely occurs because of the saturation of the CaH2 and 
CaH3 subtypes in the late-type star, which make the band- 
index classification less reliable. 

3.4. Comparison with the "Meet the Cool Neighbors" survey 

After the PMSU survey, the largest spectroscopic survey 
of M dwarfs in the northern sky was the one pr esented in the 
"Mee t the Cool Neighbo r s" (MCN) paper s e ries (ICruz & Reid l 
200 1 ICruz et alJ 12001 iReid et al.l l200l l2004t ICruz eTaLl 
2007; Reid . Cruz. & Alien! 120071) . In this section we com- 
pared the spectral clssification from the MCN survey with our 
own spectral type assignments. 

To identify the stars in common between the two surveys, 
we first performed a cross-correlation of the celestial coordi- 
nates of the stars listed in the MCN tables, to the coordinates 
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FIG. 9. — Comparison of spectral types assigned by our spectral index 
method and those assigned in the "Meet the Cool Neighbors" (MCN) spec- 
troscopic follow-up program. This shows all 161 stars in common between 
our survey and the MCN program. Spectral index subtypes are show before 
rounding up to the nearest half-integer; the MCN subtypes have random val- 
ues of ±0.2 to help in the comparison. The classifications generally agree to 
within ±0.5 subtypes (slanted lines). Our subtypes are however marginally 
later (by 0.28 subtypes on average) than the MCN subtypes. 

listed in the SUPERBLINK catalog. This was performed for 
the 1077 stars in MCN which are north of the celestial equator. 
We found counterparts in the SUPERBLINK catalog to within 
l"for 860 of the MCN stars. Of the 217 MCN stars with 
no obvious SUPERBLINK counterparts, 148 are classified as 
ultracool M dwarfs (M7-M9) or L dwarfs (L0-L7.5) which 
means they are very likely missing from the SUPERBLINK 
catalog because they are fainter than the V=19 completeness 
limit of the catalog. Of the 71 remaining stars, close exam- 
ination of Digitized Sky Survey scans failed to identify the 
stars at the locations quoted in MCN. A closer examination of 
the fields around those stars identified 49 cases where a high 
proper motion stars could be found within 3 'of the quoted 
MCN positions. These nearby high proper motion stars are 
all listed in SUPERBLINK, and have colors consistent with M 
dwarfs; we therefore assumed that the quoted MCN positions 
are in error, and matched those 49 MCN entries with the close 
by high proper motions stars from SUPERBLINK. Of the re- 
maining 22 stars we found 6 that have proper motions below 
the SUPERBLINK limit of ^ > 40 mas yr" 1 and another 5 
stars with proper motions within the SUPERBLINK limit but 
that appear to have been missed by the SUPERBLINK survey. 
Finally, there were 1 1 MCN stars that we could not identify 
at all on the Digitized Sky Survey images, and we can only 
assume that the positions quoted in MCN are too large for 
proper identification, and that the stars should be considered 
"lost". 

Of the 909 stars in the MCN program with SUPERBLINK 
counterparts, we found only 219 which satisfy the magnitude 
limit (J < 9) of our present sample of very bright M dwarfs. 
Of those, 52 stars have colors bluer (V— J < 2.7) than our sam- 
ple limit; 48 of them are classified as F and G stars in MCN, 
consistent with their bluer colors. The other 4 stars are clas- 
si fied as M dwarfs, althou gh they have V — J < 2.7 according 
to lLepine & Gaidoll d201 lh . We infer that our V -J colors for 
those stars are probably underestimated, which suggests that 
our color selection may be overlooking a small fraction of 
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nearby M dwarfs. In addition, we found another 6 stars which 
have V magnitudes and V—J colors within our survey range, 
but were rejected by t he additional infrared (J -H,H-K) 
color-color cuts used in Lepi ne & Gaidosl (1201 ll) to filter out 
red giants. All 6 stars are very bright in the infrared, and it 
appears that at least one of the H or K magnitudes listed in 
the 2MASS catalog may be in error, making the stars appear 
to have J - H and/or H-K colors more consistent with gi- 
ants. Four of the stars are classified as M dwarfs in MCN, the 
other two are late-K dwarfs. Overall, this makes a total of 8 M 
dwarfs from the MCN census that were overlooked in our se- 
lection out of the MCN subset of «150 nearby M dwarfs. This 
suggests that our color cuts, combined with magnitude mea- 
surement errors, might be missing ~ 5% of the very bright, 
nearby M dwarfs. 

In the end, this leaves only 161 stars in common between 
the MCN program and our own spectroscopic survey. The 
stars are all classified as late-K and M dwarfs by MCN, with 
subtypes ranging from K5 to M5.5. All 159 stars are identified 
with a flag ("M") in the last column of Table[3j we note that 82 
of these stars were also observed as part of the PMSU survey. 
We compare the spectral type assignments from both surveys 
in Figure|9] where the M dwarf subtypes from MCN are plot- 
ted against the (non-rounded) subtypes calculated from the 
spectral-band indices. To ease the comparison, random values 
of ±0.2 are added to the MCN subtypes. Overall, our classi- 
fications agree to within ±0.5 subtypes with the MCN values. 
The MCN subtypes, however, tend to be marginally earlier on 
average, by 0.28 subtypes; this is in contrast with the Ham- 
mer classifications (see above) which tend to be slightly later 
than our own. For the MCN subtypes, the effect is more pro- 
nounced for the earlier M dwarfs (<M2.5), where the mean 
offset is 0.43 subtypes, whereas the mean offset is only 0.09 
for the later stars. 

To investigate the difference in spectral subtype assign- 
ments, we compare the recorded values of the CaH2, CaH3, 
and Ti05 indices between the MCN program and our own sur- 
vey. After a search of the various tables published in the MCN 
series of papers, we identified 54 stars in common between 
the two programs, and for which values of CaH2, CaH3, and 
Ti05 were also recorded in both. The differences between the 
spectral index values are shown in Figure [10] For our own sur- 
vey, the corrected values of these indices are used, i.e. CaH2 c , 
CaH3 c , and Ti05 f as defined in §3.2.2. We find that the CaH2 
and Ti05 are estimated marginally higher in the MCN pro- 
gram than they are in our survey, and this very likely explains 
the difference in spectral typing: the higher index values yield 
margnially earlier spectral subtypes. This again emphasizes 
the variation in the spectral index measurements due to spec- 
tral resolution and other instrumental setups, and the need to 
apply systematic corrections between observatories to obtain 
a uniform classification system. 

We note that smaller subsets of stars in our census may 
also have spectroscopic data published in the literature, from 
various other sources. This is especially the case for the 
102 stars from the CNS3 and stars with very large proper 
motions mu > 0.2" yr" 1 which have been more routinely 
targeted for follow-up spectroscopic observations. Addi- 
tional pectroscopic survey s of selected b right M dwarfs in- 
clude^ IScholz et al] J2002l). IScholz et alJ OOOl. iReyle et al] 
(2006). and dRiaz. Gizis. & Harvinl 120051) . which all have a 
few stars in common with our catalog. Other surveys of 
nearby M dwarfs have mainly been targeting fainter stars 
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FIG. 10. — Differences between the spectral index values recorded in the 
"Meet the Cool Neighbors" (MCN) program and those measured in the 
present survey. The values are compared for a subset of 55 stars in common 
between the two surveys and for which values of CaH2, CaH3, and Ti05 were 
recorded. For the present survey, we use the corrected values (CaH2 r , CaH3 r , 
and TiO ( ) as defined in §3.2.2. The MCN values tend to be marginally higher 
on average, especially for stars of earlier spectral subtypes. This explains the 
marginally earlier spectral subtype assigments in MCN compared with the 
ones presented in this paper (see Figure|9}. 

dBochanskietalj 12501 l2010t IWest etalJl20Tlh . and do not 
overlap with our present census. 

3.5. Color/spectral-type relationships 

Spectral subt y pes were inititially estimated in 
Lepine & Gaidoi (1201 ll) based on V - J colors alone. 
Here we verify this assumption and re-evaluate the color- 
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TABLE 6 

Colors and T ef f for red dwarfs in our survey as a function of spectral 













subtype n^uv-v 1 


NUV-V 1 


cnuv-v [ 


nv-j V-J ay -J Teff 


"Teff 



K7.0 


11 


8.17 


0.21 


25 


2.90 


0.31 


4073 


98 


K7.5 


47 


8.60 


0.31 


97 


2.89 


0.16 


3883 


82 


MO.O 


87 


8.66 


0.36 


175 


2.94 


0.21 


3762 


71 


M0.5 


79 


8.74 


0.30 


159 


3.11 


0.34 


3646 


48 


M 1.0 


76 


8.89 


0.39 


151 


3.19 


0.18 


3565 


44 


M 1.5 


57 


9.07 


0.38 


147 


3.36 


0.23 


3564 


39 


M2.0 


57 


9.25 


0.47 


139 


3.52 


0.34 


3518 


57 


M2.5 


48 


9.45 


0.50 


118 


3.69 


0.28 


3500 


61 


M3.0 


34 


9.61 


0.38 


125 


3.91 


0.28 


3423 


62 


M3.5 


28 


9.69 


0.33 


124 


4.17 


0.33 


3320 


66 


M4.0 


5 


9.72 


0.35 


71 


4.45 


0.41 


3204 


76 


M4.5 


1 






36 


4.81 


0.46 


3119 


43 


M5.0 









12 


5.23 


0.50 


3014 


61 



' Non-active ("quiescent") red dwarfs only. 

magnitude relationship for bright M dwarfs. The V-J color 
index combines estimated optical (V) magnitudes from the 
SUPERBLINK catalog to the infrared J magnitudes of their 
2MASS counterparts. The SUPEBLINK V magnitudes are 
estimated either from the Tycho-2 catalog Vt magnitudes, 
or from a combination of the Palomar photographic Bj 
(IllaJ), Rf (IIIaF), and In (IVn) magnitudes, as described in 
Lepine & Shara (2005). Values of V are more accurate for 
the former (wO.lmag) than for the latter (>0.5mag); Table Q] 
indicates the source of the V magnitude. 

Mean values and dispersion about the mean of the V -J col- 
ors are listed in Table [6] for each bin of half-integer subtype; 
the table also lists how many stars of each type are in each bin. 
The V-J colors of our stars are also plotted as a function of 
spectral subtype in Figure [TTI (top panel). The ado pted color- 
subtype relationship from Lepine & Gaidosl (1201 ll) is shown 
as a thick dashed line in FigureQTJ Stars with the presumably 
more reliable Tycho-2 magnitudes are shown in blue, while 
stars with photographic V magnitudes are shown in red. Stars 
with Tycho-2 V magnitudes appear to have marginally bluer 
colors at a given subtype; this however is an effect of the vi- 
sual magnitude limit of the Tycho-2 catalog, which includes 
only the brightest stars in the V band and is t hus m ore likely 
to list bluer objects. The Lepine & Gaidos (2011) relation- 
ship generally follows the distribution at all subtypes, but with 
mean offsets up to ±0.4mag in V-J, especially at earlier and 
later subtypes. We perform a \ 2 fit to determine the following, 
improved relationship: 



[SpTy] y _ y = -32.79 + 20.75(V-J)-4.04(V-J) 2 
+0.275(V-J) 3 



(12) 



after exclusion of 3-cr outliers. The relationship is shown 
in Figure QT] (solid line). There is a scatter of 0.7 subtype 
between [SpTy] v _j and the subtype determined from spec- 
tral band indices. While the spectroscopic classification is 
more accurate and reliable, photometrically determined spec- 
tral subtypes using the equation above should still be accurate 
to ±0.5 subtype about 80% of the time, and to ±1.0 subtype 
95% of the time, which may be useful for a quick assessment 
of subtype when spectroscopic data is unavailable. 

We also compare the near-UV to optical magnitude color 
NUV — V for the 714 stars in our sample which have coun- 
terparts in GALEX; the distribution is shown in Figure QT] 
(bottom panel). We find that stars become progressively red- 
der as spectral subtype increases, from NUV - V = 8 at M0 to 



NUV-V = 10 at M4. There is however a significant fraction 
of M dwarfs which display much bluer NUV — V colors at any 
given subtype. The excess in NUV flux is strongly suggestive 
of chromospheric activity (see §6 below for a more detailed 
analysis). We separate the active stars from the more quies- 
cent objects with the following condition: 



[NUV-V] > 7.7 + 0.35 [Spty] 



(13) 



where [Spty] is the mean spectral subtype calculated from 
equations 8-11. After excluding active stars, we calculate the 
mean values and scatter about the mean of NUV-V for each 
half-integer spectral subtype. Again those are listed in Ta- 
ble [6] for reference; the table also lists the number of non- 
active stars used to calculate the mean. There is not a suffi- 
cient number of stars to calculate mean values and scatter at 
M4.5 (1 star) and M5.0 (0 star). 

4. SURVEY COMPLETENESS 

Our 1,564 spectroscopically confirmed M dwarfs are drawn 
from a catalog with proper motion limit /i < 40 mas yr -1 . The 
low proper motion limit of the SUPERBLINK catalog catches 
most of the nearby stars, but potentially overlooks nearby M 
dwarfs with small components of motion in the plane of the 
sky - either due to low space motion relative to the Sun or to 
projection effects. The catalog may also be affacted by other 
sources of incompleteness (e.g. missed detection, faulty mag- 
nitude estimate) which means that at least some very bright, 
nearby M dwarfs must be missing from our survey due to 
kinematics bias and other effects. 

To evaluate the completeness of our census, we first con- 
sider the primary source of incompleteness: the kinemat- 
ics bias of the proper motion catalog. As discussed in 
Lepine & Gaidos (201 1), the completeness depends on the lo- 
cal distribution of stellar motions, and increases with distance 
from the Sun. To estimate the kinematics bias in our sample, 
we built a model reproducing the local distribution and kine- 
matics of nearby M dwarfs. We first assumed the stars to have 
a uniform spatial distribution in the solar vicinity, and gen- 
erated a random distribution of 10 objects within a sphere 
of radius d = 70pc centered on the Sun. We assigned trans- 
verse motions to all the stars, assuming a velocity-space dis- 
tribution similar to that o f the nearby (d<10 0pc) G dwarfs 
in the Hipparcos catalog (Ivan Leeuwenll200"7l) . Because the 
distribution of stellar velocities is not isotropic, we assigned 
transverse motions for each simulated star based on the sta- 
tistical distribution of transverse motions for Hipparcos stars 
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FIG. 11. — Variation of the UV-to-optical NUV — V color index and the 
optical-to-IR V —J color index as a function of spectral subtype M. Stars with 
more reliable V magnitudes from the Tycho-2 catalog are shown in blue, stars 
with V magnitudes derived from less reliable photometric measurements are 
shown in red. The distribution of ultra-violet to optical colors (NUV— V) with 
spectral subtype shows two populations, one with a tight correlation, con- 
sistent with blackbody distribution and labeled "quiescent", and a scattered 
population of stars with clear UV excess, labeled "active". The distribution of 
V —J with subtype closely follows the relationship used by Lepine & Gaidos 



12011) to predict subtypes from V— J colors (dashed line) except for stars of 
later subtypes which have redder colors than predicted. A revised relationship 
(full line) is fitted to the data. Outliers point to stars with bad V magnitude 
measuremens. 

with sky coordinates within 30 g of the simulated object. We 
also used the simplifying assumption that the local M dwarf 
population has a uniform distribution of absolute magnitudes 
over the range 5 < Mj < 15, which is the approximate range 
of absolute magnitudes reported in the literature for M dwarfs. 

We then counted the total number of stars in the simulation 
with apparent magnitudes J < 9, and calculated the fraction of 
those stars with proper motions /i > 40 mas yr" 1 . We found 
that 93% of nearby M dwarfs with J < 9, on average, have 
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FIG. 12. — Number of spectroscopically confirmed M dwarfs as a func- 
tion of sin(b), where b is the Galactic latitude. Numbers from our census are 
shown as filled circles, with errorbars showing the Poisson noise. Distribu- 
tions predicted from models with either uniform space density, or with de- 
creases perpendicular to the Galactic plane following scale-heights of 200pc 
and 400pc, are shown for comparison. The observed distribution is largerly 
consistent with a uniform density in the local d < 70pc volume, with no evi- 
dence for incompleteness at low Galactic latitude. 

proper motions above the SUPERBLINK limit. M dwarfs 
with J < 9 extend to a maximum distance of 63 parsecs; most 
of the stars which fail the proper motion cut are in the higher 
distance range (d<50pc), and are stars near the bright end 
of the luminosity distribution Mj w 5-6. For this reason, if 
we assume a luminosity function which increases at fainter 
absolute magnitudes, the fraction of J < 9 stars which fall 
within the proper motion cut is increased, because more of 
the J < 9 stars in the local population are now M dwarfs of 
lower luminosities and closer distances, which are more likely 
to have high proper motions. The observed field M dwarf lu- 
minosity function does indeed increase for early-type dwarfs , 
to reach a peak at M r ~ 8.0 dReid. Gizis & Hawlevl I2002t 
iBochanski et al.ll20Tol) . which means that the 93% complete- 
ness estimated above must be a lower limit. To verify this, 
we tested a luminosity function where the number of stars 
increases linearly with absolute magnitude and doubles from 
Mj = 5 to My = 8; with this model, our simulations showed that 
96% of all J < 9 M dwarfs, would have /i > 40 mas yr" 1 , and 
thus be within the detection limit of SUPERBLINK. Overall, 
this suggests that only about 5% of all J < 9 M dwarfs on the 
sky will be overlooked in our census because of the proper 
motion bias. 

The SUPERBLINK surveys does however suffer from var- 
ious other sources of incompleteness, such as the inability 
to detect moving stars in saturated regions of photographic 
plates (i.e. in the immediate vicinity of very bright stars), or 
a difficulty in detecting stars in very crowded field. Also, the 
SUPERBLINK code has trouble detecting the motions of rel- 
atively bright V < 12 stars because of the saturated cores of 
their point spread functions. In practice, this is mitigated by 
incorporating data from the Tycho-2 catalog, which provides 
very accurate proper motion measurements for bright stars. 
However the Tycho-2 catalog itself has some level of incom- 
pleteness in the 8 < V < 12 manitude range. 
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One way to test incompleteness due to crowding or satu- 
ration is to examine the distribution of SUPERBLINK stars 
as a function of Galactic latitude. Crowding and saturation 
effects should be more pronounced in low Galactic-latitude 
fields, where the stellar density is high. In comparison, stars 
in high Galactic latitude fields will be easier to detect with 
the SUPERBLINK code. The same is true for the Tycho-2 
catalog, which should be more complete at high Galactic lati- 
tudes. 

We calculated the number of spectroscopically confirmed 
M dwarfs in our sample as a function of sin(b), where b is the 
Galactic latitude. The distribution is shown in Figure [TU the 
number of stars in each bin is plotted as a filled circle, with 
errobars showing the 1 a Poisson error. The distribution in- 
creases with sin(b) because our stars are all located north of 
the celestial equator. For comparison, we plot the trend ex- 
pected of a uniform distribution of stars in the local volume 
(blue histogram), assuming the same total number of stars as 
in our census. We find our data to be largely consistent with 
the uniform distribution model, and see no evidence of a sig- 
nificant dip at low Galactic latitude, which one would expect 
if the SUPERBLINK survey is incomplete at low b. This sug- 
gests that our sample does not suffer from significant sources 
of incompleteness due to saturation/crowding. 

One might however argue that the uniform density assump- 
tion is invalid for a census extending to ~70pc of the Sun, 
and that a slight overdensity of stars should be detected at low 
Galactic latitudes. The fact that we see no evidence of such 
an overdensity in our data would then be indirect evidence 
of incompleteness at low Galactic latitude. To examine this 
possibility, we estimated the expected distributions from sim- 
ulations in which the stellar density perpendicular to the plane 
(i.e. along Z) decreases with scale heights of 400pc or 200pc. 
Kinematics and absolute magnitude distributions were also 
simulated as described above, and stars were selected based 
on J < 9 and /i > 40 mas yr" 1 . The mean distributions from 
the 400pc and 200pc scale-height models are shown in Fig- 
ure [12] and do indeed predict a slight overdensity of objects 
at low Galactic latitudes. Assuming that all the stars at high 
Galactic latitude are detected in all the models, we find a ^5% 
excess of stars in the 400pc scale-height model over the ob- 
served number, and a ~10% excess of stars in the 200pc scale- 
height model. One might then argue that the SUPEBLINK 
census potentially has an additional 5%-10% incompleteness 
level, due to incompleteness in the plane of the Milky way. 
This is most probably an overestimate, however, because there 
is no evidence that the local stellar distribution shows any sig- 
nificant decrease with Z. Overall, the SUPERBLINK census 
appears to be essentially complete at low Galactic latitudes. 

Proper-motion selection may introduce an additional bias 
against metal-rich and/or young stars, which tend to have 
lower components of motion in the vicinity of the Sun. How- 
ever this effect is expected to be small, e.g. a few per- 
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(Perrvman et al.1119981). U rsa Majoris (King et al.ll2003l) . and 
TW Hvdrae dReidl 120031) young nearby moving groups have 
proper motions that exceed 40 mas yr" 1 and thus would not 
be selected against using the current selection method. 

Finally, some bright M dwarfs may be overlooked because 
of the imposed color cuts, due to magnitude errors and uncer- 
tainties. As described in §3.4 above, a handful of previously 
known M dwarfs from the MCN census were overlooked in 
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FIG. 13. — Effective temperatures for the bright M dwarfs in our survey 
determined by fits to the PHOENIX models. The gray points represent in- 
dividual objects, while the large black points are the median values of the 
objects within each half subtype bin. The error bars are the interquartile 
ranges. Blue diamonds and red triangles show the T^ff estimates for sub- 
sets of stars in our survey whose tempreatures were estimated from photom- 
etry Casagrande, Flynn, & Bessell (2008) and model fits to infrared spectra 
iRoias-Ayala et al. 2012) respectively. While both spectroscopic estimates 
show evidence of a plateau at M2-M3, the photometric estimates do not con- 
cur. 

our target selection for precisely those reasons. As suggested 
in §3.4, it is possible that we may be overlooking 5% of bright 
M dwarfs because of magnitude errors. This, combined with 
the estimated 96% completeness in the proper motion selec- 
tion estimated above, suggests that our list of 1405 M dwarfs 
likely includes w 91% of all existing M dwarfs with infrared 
magnitude J < 9 as seen from Earth. There may still be rs 140 
bright M dwarfs to be identified, although the precise number 
can only be determined after these "missing" stars are found. 

5. PHOENIX MODEL FITS AND T EFF ESTIMATES 

We compared our spectra to a grid of 298 models of 
K- and M-dwarf sp ectra genera ted by the BT-SETTL ver- 
sion of PHOENIX dAllard et al.ll2010l) . BT-SETTL includes 
updated opacities (i .e. of H2O), revised solar abundances 
( Aspl und et al.120091) . a refractory cloud model, and rotational 
hydrodynamic mixing. The models include effective tempera- 
tures T eff of 3000-5000 K in steps of 100 K, logg values of 4, 
4.5, and 5, and metallicities of [M/H] = -1.5, -1, -0.5, 0, +0.3, 
and +0.5. For each temperature, logg and metallicity value, 
we selected the model with a/Fe that was closest to solar. 

The spectral density of model calculations varies with 
wavelength but is everywhere vastly greater than the resolu- 
tion of our spectra. Model spectra were thus convolved with 
a gaussian with FHWM of the same resolution as the spectra, 
and a corrective shift (typically less than a resolution element) 
was found by cross-correlating the observed and model spec- 
tra. Normalized spectra were ratioed and \ 2 calculated us- 
ing the variance spectrum of the observations. We restricted 
the spectral range over which \ 2 is calculated to 5600-9000A 
and excluded the problematic region 6400-6600A which con- 
tains poorly-modeled TiO absorption dReyle et al.ll201 ll) . We 
also excluded regions where the telluric correction is rapidly 
changing with wavelength, i.e. the slope, smoothed over 4 
resolution elements or 11.7 A, exceeds 1.37 A" 1 . The model 
with the smallest \ 2 value was identified. For a more refined 
estimate of effective temperature, we selected the 7 best-fit 
models and constructed 10,000 linear combinations of them; 
the "effective temperature" of each is the weighted sum of the 
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FIG. 14. — Distribution of corrected CaH2+CaH3 vs Ti05 band index values for the stars in our survey (grey dots - with brightness levels correlating with 
spectral subtype). There is a tight correlation between the two indices, which are both also correlated with spectral subtypes, with earlier stars on the upper 
right of the diagram as shown. The distribution is used as a guide to calibrate the value of £, with C = 1 assu med to trace the CaH2+CaH3/T i05 relationship 
for sta rs with average Galactic disk abundances (near-solar). The iso-f contours from the earlier calibrations of Lepine, Rich, & Shara (2007) and Dhital et al. 
1 2012) are shown as dotted a nd dashed lines, respec t ively . When applied to our corrected spectral index values, they both d iverge from the obs erved distribution 
at earlier subtypes, with the Lepine, Rich, & Shara (2007) overestimating the £ of late-K and early-M dwarfs, while the D hital et al. 12012) calibration yield 
underestimates. This emphasizes again the need to use properly recalibrated and corrected spectral index values (see Figures l4l5t . Our revised, dataset-specific 
calibration is shown with the continuous lines. Known metal-rich amd metal-poor stars are denoted in red and blue, respectively. 

temperatures of the components. We again found the model 
with the minimum \ 2 - We calculated the standard deviation of 
T e ff among the combination models as a function of the max- 
imum allowed x 2 - We reported the maximum standard devia- 
tion as a conservative estimate of uncertainty. We also calcu- 
lated formal 95% confidence intervals for T e ff based solely on 
the expected distribution of \ 2 f° r N — 3 degrees of freedom, 
where N ~ 1 100 is the number of resolution elements used in 
the fit. The parameters of the best-fit model, and the refined 
T e ff, standard deviation, and confidence intervals are reported 
in Tabled 



Values of T e /f calculated for individual stars are plotted in 
Figure [13] as a function of their spectral subtype (grey dots) . 
Median values for stars within each half-subtype bin are plot- 
ted in black, with error bars showing the interquartile ranges. 
Our model-fit algorithm prefers values that match the T e ff 
model grid, which have a 100K grid step (i.e., 3500K is pre- 
ferred over 351 OK). 

Our results suggest the existence of a T e ff plateau 
spanning M1-M3. To investigate this further, we com- 
pare our values to the effective temperatures reported in 
Casagrand e. Flvnn. & Besselil (2008) for 18 of the stars in our 
sample, and in dRojas-Ayala et alJl2012l) for another 49 stars ; 
our own spectral type determinations are combined to the T e ff 
measured by the other authors. The values are compared in 
Figure [13] We find that our mid-type plateau is corroborated 



with the (Roias-Avala et al. 2012) values, but not with those 
fromlCasagrand e. Flvnn. & Bessell (2008), whose values de- 
crease more linearly with spectral subtype. The eff ective tem- 
peratures in Casagrande, Flynn, & Bessell (200 8]) are based 
on ph otometric measurements while the (Rojas-Ayala et al. 
120121) values are estimated by PHOENIX model fits to in- 
frared spectra. It is interesting that the fits to the optical and 
infrared spectra yield T e ff values which are in agreement. The 
disagreement with the photometric determinations however 
suggest that atmospheric models for M dwarfs are still not 
well understood, in particular in the M1-M4 spectral regimes. 

6. THE f -PARAMETER AND METALLICITY ESTIMATES 

6.1. Recalibration of the ( parameter 

The Cno/CaH parameter (denoted ( for short) is a combina- 
tion of the Ti05, CaH2, and CaH3 spectral indices which was 
shown t o be correlated with metallicity in m etal-poor M sub- 
dwarfs (IWoolf. Lepine. & Wall erstein 2009). The index was 
first described in Lepine, Rich, & Shara (200 7|), and a revised 
calibration has recently been proposed by lDhital et al.l(l2012l) . 
The index measures the relative strength of the TiO molecular 
band around 7,000A with respect to the nearby CaH molecu- 
lar band. In cool stars, in fact, the ratio between TiO and CaH 
is a function of both gravity and metallicity. The CaH band 
is noticeably stronger in giants dMann et al.ll2012l) . and this 
effect can be used as affective means to separate out M giants 
from M dwarfs using optical spectroscopy. In the higher grav- 
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FIG. 15. — The £ parameter as a function of spectral subtype. Top: differ- 
ence in f for the same stars measured at the two observatories (MDM and 
UH). The scatter provides an estimate of the measurement error on ( , which 
is significantly larger at earlier subtypes. Bottom: adopted values of ( for all 
the stars in the survey. The larger scatter at subtypes M2 and earlier can be 
fully accounted by the measurement errors. The scatter at subtypes M3 and 
later is larger than the measurement error, and is thus probably intrinsic and 
is evidence for intrinsic metallicity scatter in the solar neighborhood. 

ity M dwarfs/subwarfs however, the TiO to CaH ratio is how- 
ever believed to be mostly affected by metallicity, although 
young stars may show gravity effects as well. 

The £ parameter was originally introduced to rank metal- 
poor, main-sequence M star s into three metallicity classes 
dLepine. Rich, & Sharall2007l) : stars with 0.5 < C < 0.825 are 
formally classified as subdwarfs (sdM), 0.2 < ( < 0.5 de- 
fines extreme subdwarfs (esdM), while a ( < 0.2 identifies 
the star as an ultrasubdwarf (usdM). However, it is conjec- 
tured that £ could be used to measure metallicity differences 
in disk M dwarfs, i.e. at the metal-rich end. Disk M dwarfs 
are generally found to have 0.9 < ( < 1.1, though it is un- 
clear if variations in £ correlate with metallicity for values 



within that range. Measurement of Fe lines in a subset of 
M dwarfs and subdwarfs does confirm that the C parameter 
is cor related with metallicity (IWoolf. Lepine. & Wallersteinl 
2009), with £ — 1 .05 presumably corresponding to solar abun- 
dances. However there is a significant scatter in the relation- 
ship which raises doubts about the accuracy of £ as a metal- 
licity diagnostic tool. 

A important caveat is that the TiO/CaH ratio is not sen- 
sitive to the classical iron-to-hydrogen ratio Fe/H, but rather 
depends on the relative abundance of a-elements to hydro- 
gen (a/H) because O, Ca, and Ti are all a-elements. Vari- 
ations in a/Fe would thus weaken the correlation between 
£ and Fe/H. The a/Fe abundance ratio is however relatively 
small in metal-rich stars of the thin disk (±0.05dex) disk 
stars, and are found to be significant (±0.2dex) mostly in 
more metal-poor star s associated with the thick disk and halo 
dNavarro et al.ll2.0lTb . It is thus unclear whether typical a/Fe 
variations would affect the £ parameter significantly in our 
subset, which is dominated by relatively metal-rich stars. 

On the other hand, it is clear that the index has significantly 
more leverage at later subtypes. This is because the strengths 
of both the TiO and CaH bands are generally greater, and their 
ratio can thus be measured with higher accuracy. The index 
is much leass reliable at earlier M subtypes however, and is 
notably inefficient for late-K stars. 

A more important issue is the specific calibration adopted 
for the £ parameter, which is a complicated function of the 
Ti05, CaH2, and CaH3 indices. The £ parameter itself is 
defined as: 

^-Ti05_ 
^ l-[TiO5] Z0 ' 

which in turns depend on [TiO5]z , itself a function of 
CaH2+CaH3. The function [TiO5]z represents the expected 
value of the Ti05 index in stars o f solar metallicity, for 
a giv en value of CaH2+CaH3. In (Lepi ne. Rich. & Sharal 
120071) . Ti05] Zo was defined as: 

[Ti05] Zp = -0.050 - 0. 1 1 8 CaH+ 0.670 CaH 2 



-0.164 CaH 



(15) 



where CaH = CaH 2 + CaH3. The more recent calibration of 
Dhit al et ail (120121) . on the other hand, uses: 

[TiO5]z = -0.047-0.127 CaH+0.694 CaH 2 

-0.183 CaH 3 -0.005 CaH 4 . (16) 

The difference between the two calibrations is mainly in the 
treatment of late-K and early-type M dwarfs, as illustrated in 
Figure [14] When overlaid on the distribution of CaH2+CaH3 
and Ti05 values from our current survey, however, the two 
calibrations fail to properly fit the distribution of data points at 
the earliest subtypes (hig h values of CaH2+CaH 3 and Ti05). 
This results in the Lepine, Rich, & Shara (20071) overestimat- 
ing the metallicity at earlier subtypes, while the iDhital et"aT] 
(2012) calibration tends to underestimate metallicity. 

In any case, the evidence presented in §3.3 and §3.4 which 
shows that differences in spectral resolution and flux cali- 
bration can yield differences in the Ti05, CaH2, and CaH3 
spectral indices of the same stars, also suggests that a cali- 
bration of the £ parameter may only be valid for data from 
a particular observatory /instrument. A general calibration of 
£ may only be adopted after corrections have been applied 
as described in Section 3.2. Because we do not have any 
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star in common with the iDhital et al.l (120121) subsample, we 
cannot verify the consistency of their £ calibration to our 
data at this time. In addition, because we have now applied 
a correction to our MDM sp ectral index measurements, the 
iLepine. Rich. & Sharal d2007) calibration of £ may now be 
off, and should not be used for our sample. 

Instead, we recalibrate the £ parameter again, using our cor- 
rected spectral index values. Our fit of [Ti05] c as a function 
of [CaH] c = [CaH2] c + [CaH3] e yields: 

[TiO5]z =0.622- 1.906 ([CaH] c )- 2.211 ([CaH] c ) 2 

-0.588 ([CaH] c ) 3 . (17) 

We calculate the new £ values using the corrected values of 
the Ti05 index, i.e.: 



l-[Ti05] c 
l-[Ti05] Z(T 



(18) 



All our values of £ are listed in Table [3] 

In order to evaluate the accuracy of the £ measurements, 
we compared the values of £ independently measured at both 
MDM and UH for the 146 stars in our inter-observatory sub- 
set. Values are compared in Figure Q3] (top panel) which 
shows A£ = (mdm - Qjh as a function of spectral subtype. We 
find a mean offset A£ = -0.01 and a dispersion ctac = 0.10. 
The small offset indicates that the £ measurements are gener- 
ally consistent between the two observatories. The dispersion 
<ta(, on the other hand, provides an estimate of the measure- 
ment accuracy. Splitting the stars in three groups, we find the 
mean offsets and dispersions (A£,<tac) t0 ^ e (-0.086,0.244) 
for subtypes K7.0-M0.5, (-0.011,0.103) for subtypes M1.0- 
M2.5, and (0.008,0.036) for subtypes M3.0-M5.5. Assuming 
that stars do not show significant variability in those bands, 
we adopt the dispersions as estimates of the measurement er- 
rors on £ for that particular subtype range. It is clear from 
Figure [14] that early type stars should have larger uncertain- 
ties in £ because of the convergence of the iso-£ lines. The 
best leverage for estimating metallicities from the TiO and 
CaH bandheads is at later types when the molecular bands 
and well developed. 

The overall distribution of £ values as a function of spectral 
subtype also shows a decrease in the dispersion as a func- 
tion of spectral type (Figure [15] bottom panel). In early-type 
dwarfs (K7.0-M0.5), the scatter in the £ values is relatively 
large, with ~ 0.174. It then drops to ~ 0.100 for sub- 
types Ml .0-M2.5, and to ct c ~ 0.059 for subtypes M3.0-M5.5. 
Note that the scatter in the M3.0-M5.5 bin is a factor 2 larger 
than the estimated accuracy of the £ for that range, as esti- 
mated above. We suggests this to be evidence of an intrin- 
sic scatter in the £ values for the stars in our sample, which 
we allege to be the signature of a metallicity scatter. If we 
subtract in quadrature the 0.035 measurement error on £, we 
estimate the instrinsic scatter to be «0.05 units in £. This in- 
trinsic scatter, which presumably affects all subtypes equally, 
is unfortunately drowned in the measurement error at earlier 
subtypes (<M3). 

6.2. Comparison with other metallicity diagnostics 

To test our £ as a tracer of metallicity for dM stars from M3 
to M6, we compared the values to two recent [Fe/H] calibra- 
tion techniques for M dwarfs with solar metallicities . First , 
we used the photometric calibration of [Neves et ail (120121) . 
which is based on optical-to-infrared V-K color and absolute 



0.6 
0.4 
2 0.2 
5 0.0 
^ -0.2 
-0.4 
-0.6 





0.6 




0.4 


CM 
< 

cc 


0.2 




0.0 








-0.2 




-0.4 




-0.6 



0.9 



M3 
M4 
M5 
M6 



1.0 1.1 



1.2 



0.9 



1.0 



1.1 



1.2 



FIG. 16. — Comparison between the £ parameter values and independent 
metallicity measurements for subsets of M dwarfs in our survey. Top panel 
compares our ( to the Fe/H estimated from the (V-K.M^) calibration of 
INeves et al. 1 2012) for the same stars. Bottom panel com pares our C values 
to the Fe/H estimated from the infrared K-band index by Roias-Avala et al. 
1 2012). Both distribution show weak positive correlations. 

magnitude M^. The method is sensitive to small variations in 
V-K/Mk and thus requires an accurate, geometric parallax. A 
total of 143 stars in our sample have parallaxes, and thus can 
have their metallicities estimated with the method. Figure [161 
(top panel) plots the estimated [Fe/H] as a function of £ for 
those 143 stars. The distribution shows significant scatter, but 
we find a weak correlation of [Fe /H] with £, which we fit with 
the relationship: 



[ J Fe/i/] N12 = 0.750C-0.743. 



(19) 



Stars are scattered about this relationship with a I- a disper- 
sion of 0.383 dex. One drawback of the photometric metallic- 
ity determination is that it assumes the star to be single. Unre- 
solved double stars appear overluminous at a given color, and 
will thus be determined to be metal-rich. Also, young and ac- 
tive stars often appear overlu minous in the color-magnitude 
diagram dHawlev et al.ll2002h . and their metalicities based on 
V-K/M/s- would also be overestimated. Multiplicity and ac- 
tivity could therefore contribute in the observed scatter. Stars 
with [Fe/H]^i2 > 0.4, in particular, could be overluminous in 
the V-K/Mk diagram, as their £ does not suggest them to be 
metal-rich. 

Next, we re tr ieve metallicity measurements from 
iRoias-Avala et al.l d2012l) . who estimated [Fe/H] based 
on the spectroscopic calibration from infrared K-band atomic 
features. Their list has 37 stars in common with our survey. 



Spectroscopic catalog of the brightest, northern M dwarfs 



21 



:0.60 0.70 0.80 0.90 


i i i i i 
00 1 10 


i 1 




1 .0 


^^^^^^ 


1 


. • • 




0.8 


f values 




• • • . • : 

■ 

/; 




0.6 


• * •"' v . 

* " * 




0.4 


• » 






_D 


0.2 


... 

• * "V. 

■ . • • V. tdK 

. " • *! \j**Pff\ 






~a> 

Q 


0.0 














i 






-0.2 
-0.4 



-2 2 4 

sptype 




2 
sptype 



2.0 




2.0 



3.0 



.6 



A 



.0 
0.5 




BT-SETTL log g= 5 
1.0 1.5 2.0 



2.5 



3.0 



FIG. 17. — SDSS photometry of M dwarf stars synthesized from SNIFS spectra and transmission functions convolved with unit airmass. Upper left: r-Ks 
(where Ks is from the 2MASS point source catalog) against M spectral subtype (where K7=-l and K5=-2). Points are colored by the £ parameter, which 
measures TiO/CaH ratio and is a metallicity diagnostic in the optical. The f values are undefined for late K stars, which are plotted as black points. The black 
curve is a running median (N=81). Upper right: difference of r-Ks with respect to the running median vs. spectral type, showing no obvious co rrelation w ith zeta. 
Lower left: g-r vs. r-z, showing an apparent correlation between these colors and £. The contours are the empirical function f or £ derived by West et al .1 )201 II) . 
Lower right: SDSS g-r vs. r-z colors generated by the PHOENIX/BT-SETTL model I Allard, Hauschildt, & Schweitzer 200jj) for log g=5, Teff=3500K-4200K, 
and five different values of the metallicity as noted in the legend. The model predicts the more metal-rich M dwarfs to have bluer g-r colors, while being redder 
in r-z. The color dependence on metallicity in most pronounced in late-type stars, and nearly vanishes at K7/M0. 



The [Fe/H] values are plotted as a function of our £ values 
in the bottom panel of Figure [16] Again there is significant 
scatter, but we also find a weak correlation which we fit with 
the relationship: 



[Fe/H] RA12 =l.07l(- 1.096, 



(20) 



about which there is a dispersion of 0.654 dex. The statistics 
are relativey poor at this time, and more metallicity measure- 
ments in the infrared bands would be useful. 

The weak correlation found in both distribution is interest- 
ing in itself. Using a sample of stars spanning a wide range of 
metallicities and £ values, including metal-poor M subdwarfs, 
extreme subdwarfs (esdM), and extremely metal-poor u ltra- 
subdwarfs (usdM), IWoolf. Lepine. & Wallersteinl (120091) de- 
termined a relationship of the form [Fe/H] = -1.685+ 1.632£, 
over the range 0.05 < C < 1-10. All the stars in the two distri- 
butions from the present survey have ( values between ^0.9 



and ~1.2, and thus represent the metal-rich end of the distri- 
bution. The weaker slopes we find in our correlations (0.75 
and 1.07) may indicate that the relationship levels off at high 
metallicity end, which would make £ much less useful as a 
metallicity diagnostic tool in Solar-metallicity and metal-rich 
M dwarfs. The correlations are however weak, and more ac- 
curate measurement of £ and Fe/H would be needed to verify 
this conjecture. 

6.3. The one M subdwarf: PM 120050+5426 (VI 51 3 Cyg) 

The primary purpose of the £ parameter is the identifica- 
tion of metal-poor M subdwarfs, for which it has already 
proven ef fective. By definition, M sub dwarfs are stars with 
C <0.82 dLepine. Rich. & Sharal 120071) . Though we have a 
few stars with values of £ just marginally under 0.82, only 
one star clearly stands out as a definite M subdwarf: the star 
PM 120050+5426 (= Gl 781) which boasts a C = 0.58 well 
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within the M subdwarf regime. The star also clearly stands 
out in Figure[l4]where it lies noticeably below the main locus 
at Ti05 c -0.75. 

PM 120050+5426 is also known as V1513 Cyg, a star pre- 
viously identified as an M subdwarf by Gizis ( 1997), a nd one 
clearl y associated with the Galactic halo dFuchs & Jahreissl 
1998). The star is also notorious for being a flare star, with 
chromospheric activity due not to young age bu t to the pres - 
ence of a low-mass companion on a close orbit (Gizis 1998). 
Our own spectrum indeed shows a relatively strong line of Ha 
in emission, which is extremely unusual for an M subdwarf. 
It is an interesting coincidence that the brightest M subdwarf 
in the northern sky should turn out to be a peculiar object. 

In any case, because the TiO molecular bands are weaker 
in M subdwarfs than they are in M dwarfs, the use of TiO 
spectral indices for spectral classification leads to underesti- 
mates of their spectral subtype. The convention for M sub- 
dwarfs is rather to base the clas sification on the strengths of 
the C aH bandheads dGizisll 19971: ILepine. Rich. & Sharal200"l 
120071) . We adopt the same convention here, and recalculate 
the subtype from the mean of Equations 6 and 7 only (CaH2 
and CaH3 indices). We thus classify PM 120050+5426 as an 
sdM2.0, which is one h alf-subtype l ater than the sdM1.5 clas- 
sification suggested by (Gizis 1997). 

6.4. Photometric dependence on metallicity 

A prediction of current atmospheric models is that metallic- 
ity variations in M dw arfs yield significa nt variations in opti- 
cal broadband colors (lAllard et al.ll2000l) . The metal-poor M 
subdwarfs have in fact long been known to have bluer V-I col- 
ors tha n the more metal-rich field M dwarfs of the same lumi- 
nosity dMonet et al.llT99l iLepine. Rich. & Sharall2003l) . The 
bluer colors are due to reduced TiO opacities in the optical, 
which make the spectral energy distribution of M subdwarfs 
closer to that of a blackbody, while it makes the metal-rich M 
dwarfs display extreme red colors. 

Interestingly, the SDSS g — r color index shows the oppo- 
site trend, and is bluer in the more metal-poor stars. This is 
because the TiO bands very strongly depress the flux in the 
6000A-7000A (r-band) range, an effect which in fact makes 
the metal-rich M dwarfs degenerate in g-r, as the increased 
TiO opacities in cooler stars balance out the reduced flux in g 
from lower T e yy. This effect is much weaker in metal-poor 
stars due to the reduced TiO opacity, which makes metal- 
poor stars go redder as they are cooler, as one would nor- 
mally expect. This has been observed in late-type M sub- 
dwarfs, which have sign i ficant ly redder color that field M 
dwarfs ( ILepine & Scholzl 120081) . The color dependence of 
M dwarfs/subdwarfs on metallicity is also predicted by at- 
mospheric models. Figure [T7l( bottom-right panel) shows the 
predicted g-r and r-z colors from the PHOENIX/ BT-SETTL 
model of Allard, Hauschi ldt. & Schweitzerl (12000). The mod- 
els corroborate observations and predict redder g-r colors in 
metal-poor stars. 

Although ugriz photometry is not available for our stars 
(all of them are too bright and saturated in the Sloan Digital 
Sky Survey), it is possible to use the well-calibrated SNIFS 
spectrophotometry to calculate synthetic broadband riz mag- 
nitudes for the subset of stars observed at UH. We first ex- 
amine any possible correlation between the optical to infrared 
r-Ks color (taken as a proxy for V-I) and the £ index. Fig- 
ure[l7]plots r-K s as a function of spectral subtype, with the 
dots color-coded for the ( values of their associated M dwarf 
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FIG. 18. — Fraction of active stars as a function of the spectral subtype M. 
The rise at later subtypes is consistent with earlier studies of field M dwarf, 
which shows increased activity levels in mid-type M dwarfs. The fraction 
level at later subtype is however higher than that measured in the SDSS spec- 
troscopic catalog. 

(top-left panel). We find a tight relationship between r-K s 
and spectral subtype, which we fit using a running median. 
The residuals are plotted in the top-right panel, and show no 
evidence of a correlation with (. There are a significant num- 
ber of outliers with redder r-K s colors than the bulk of the M 
dwarfs: these likely indicate systematic errors in estimating 
the synthetic r band magnitudes. The absence of any clear 
correlation suggests that an optical-to-infrared color such as 
r-Ks is not sensitive enough to detect small metallicity vari- 
ations, at least at the metal-rich end. 

The synthetic g-r and r-z colors are plotted in Figure HTl 
(lower-left panel). The redder stars (r-z>1.2) show a wide 
scatter in g-r, on the order of what is predicted for stars with 
a range of metallicities -0.5 < [Fe/H] < 0.5. Though we do 
not find a clear trend between the synthetic g-r colors and 
the £ values measured in the same stars, the high-C stars (red 
and orange dots on the plot) do seem to have lower values of 
g-r on average than the low-C ones (green dots). The trend 
is suggestive of a metallicity link to both the g-r colors and 
the C values, and should be investigated further with data of 
higher precision. 

7. CHROMOSPHERIC ACTIVITY 

To evaluate the presence of Ha in emission in our M dwarfs, 
we used the Ha equivalent width index EWHA defined as: 



EWHA = 100 



100 



0(A)rfA+0(A>*A 



(21) 

where S(\) is the observed spectrum. The EWHA index mea- 
sures the flux in a region (6557.61A-6571.61A), which in- 
cludes the Ha line, in relation to a pseudo-continuum region 
spanning 6500A-6550A and 6575A-6625A; the calculation 
provides a value in units of wavelength (A) like the traditional 
equivalent width. Note that for an Ha line in emission, values 
of the EWHA index are negative, following convention. As- 
suming the W1-W2 region to measure the true spectral con- 
tinuum, then the EWHA index would measure the true equiv- 
alent width of Ha. As it turns out, the Wl - W2 region often 
has a higher mean flux than the W3 — W4 region without the 
Ha emission component, which means that the EWHA index 
systematically underestimate the equivalent width of the Ha 
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TABLE 7 

Survey stars: distances, kinematics, and activity. 



Star name 



V>or 



U 
km s _1 



V 
km s 



W 
km s 



EWHA 
A 



Ha 
active 



Xray 
active 



UV 
active 



100006+1829 

I00012+1358S 

100033+0441 

100051+4547 

100051+7406 

100077+6022 

100078+6736 

100081+4757 

100084+1725 

100088+2050 

100110+0512 

100113+5837 

100118+2259 

I00125+2142En 

100131+7023 



0.0342± 0.0032 
0.0889± 0.0014 



0.0460± 0.0019 



0.0233± 0.0038 



0.0358± 0.0028 



0.030± 0.008 


0.028± 0.008 


-13.9 


12.1 




0.41 






0.031 ±0.008 


0.030± 0.009 


8.5 


-6.8 




0.39 






0.078± 0.021 


0.083± 0.024 


-38.2 




-15.9 


0.47 






0.078± 0.031 


0.091 ±0.027 


-15.1 




-4.4 


-3.11 


Y 




0.046± 0.012 


0.055± 0.016 


5.1 




-8.5 


0.39 






0.071 ±0.028 


0.061± 0.018 


7.9 




1.8 


-2.93 


Y 


Y 


0.040± 0.01 1 


0.040± 0.012 


11.2 


0.7 




0.41 






0.078± 0.031 


0.080± 0.024 


9.5 




-10.1 


-4.64 


Y 


Y 


0.032± 0.009 


0.031 ±0.009 


-48.5 


-14.3 




0.35 






0.055± 0.022 


0.054± 0.016 


-3.0 




-16.7 


0.30 




Y 


0.023± 0.006 


0.024± 0.007 


-5.9 




-32.5 


0.44 






0.037± 0.010 


0.037± 0.011 


-6.2 




16.8 


0.37 







PM 
PM 
PM 
PM 
PM 
PM 
PM 
PM 
PM 
PM 
PM 
PM 
PM 
PM 
PM 



line. The index is however reproducible and more convenient 
for automated measurement than, e.g. manual evaluation of 
the equivalent using interactive software such as IRAF. 

The EWHA index was measured for all spectra in our sam - 
ple, and used to flag active stars. Following West et al.l (1201 ll) . 
we defined a star to be chromospherically "active" if EWHA 
< -0.75, which usually corresponds to a clearly detectable 
Ha line in emission. Values of the EWHA index are listed 
in Table [7] Under the above criterion, 171 M dwarfs in our 
s urvey are considered act i ve. 

lHawlev. Gizis. & Reidl d 19961) found that active stars (by 
their criterion, EW(Ha) > 1) have slightly redder V -K col- 
ors for the same value of Ti05 index, an effective tempera- 
ture proxy. We calculated a median (N = 30) locus of V - K 
vs. Ti05 index for our entire sample and found that 13 active 
stars are bluer than this locus, while 45 are redder, seemingly 
confirming their result. The large scatter in V - K colors how- 
ever prevents us from quantifying this offset more precisely. 

The fraction of stars that are active at each spectral subtype 
is shown in Figure [18] with error bars computed from the bi- 
nomial distribution. The increase in the active fr action with 
spect r al type is consistent with previous studie s dJov & Abtl 
19741: lHawlev. Gizis. & Reidl 1 19961. IWest et alJ 120041 120081: 
Kruse et all 120101: IWest et alJ 1201 ll) . Our active fractions 
are hig her at subtype M4-M6 than the lHawlev. Gizis. & Reidl 
( 1996) results, even when using their criterion to define active 
stars (see above). This may be a result of a slightly differ- 
ent definition for EW, or a result of the relatively small num- 
ber of objects. O ur active fractions are also higher than the 
IWest et all (1201 ll) results at each subtype. This discrepancy is 
likely caused by the magnitude limit imposed in our survey: 
our objects are all nearby, and relatively close to the Galactic 
plane (see Sectio n 8), which makes t hem statistically younger, 
as also suggested IWest et ail (1201 ll) . The active f ractions for 
our sta rs are closer to the active fractions for the IWest et aP 
(1201 ll) stars in bins of stars closest to the Galactic plane. 

Another chromospheric activity diagnostic in M dwarfs 
is the detection of X-rays. M dwarfs that are X-ray 
bright are often young, and this has been used t o iden- 
tify members of nearby young moving groups (e .g., [Gaidos 
119981: IZuckerman et all 12001: Tor res et al.l 120061). and other 
young stars in the Solar N eighborhood (|Riaz. Gizis. & Harviiil 
2005). Most recent ly, IShkolnik et al.l (120091 120121) and 
Schlieder etakl (120 12b have used the ratio of ROSAT X-ray 
flux to 2MASS J or /T-band flux to identify candidate mem- 



bers of young moving groups. This technique is particu- 
larly effective for objects < 70pc away, which inc ludes all 
the stars in our sample. The Lepine & Gaido| (1201 ll) catalog, 
from which our targets are drawn, was already cro ss-matched 
to the ROSAT All-Sky Bright Source Catalog (Vog es et al l 
119991) and the RO SAT All-Sky Survey Faint Source Catalog 
(Vog es et alfcOOOl) . 

We have computed the X-ray flux for our survey stars from 
the measured c ount rate and hardness ratio (HR1) us ing the 
prescription in Schmitt, Fleming, & Giampapa (1995). Fig- 
ure [19] shows the distribution of X-ray flux as a function of 
V-K color, for the 290 M dwarfs with ROSAT detections. 
Dots are color-coded according to the strength of the Ha 
emission, as measured by the EWHA index. As expected, 
M dwarfs with strong Ha emission also tend to be more X- 
ray bright. Objects with log Fx Fk > -2.6 (above the dashed 
line in Figure [i9l are considered bright enough in X-ray to 
qu alify as chromo s pheric ally active, following the definition 
of Schliederitin (I20T2I) . Some 154 of the X-ray sources 
are active based on their EWHA values, and all of them also 
qualify as active stars based on their X-ray fluxes. On the 
other hand, 53 M dwarfs identified as active based on X-ray 
flux do not display significant Ha emission in our spectra; 
most of them tend to be earlier M dwarfs, in which Ha emis- 
sion is not as easily detected as in later type objects because 
of their higher continuum flux near A6563. There are also 
22 stars in our survey which are active based on Ha but are 
not detected by ROSAT. This suggests that only two thirds 
of the "active" stars will be diagnosed as such from both X- 
ray and Ha emission, while the other third will show only 
either. This could be due to source confusion in the ROSAT 
X-ray survey, variability in either X-ray or Ha emission, or, 
in the case of the X-ray flux, non-uniform sky coverage by 
ROSAT. We calculated the luminosity ratio index Lx /Ln a of 
active stars, as d efined by the criterion EW(Ha) > 1 following 
the procedure of lHawlev. Gizis. & Reidl (119961) . and adopting 
the relation V -R xs 0.7 + 0.06SpTy to estimate an R mag- 
nitude and the continuum flux at the Ha line. We find that 
the ratio is insensitive to bolometric magnitude and spectral 
type, and has a median value of 0.85. This is higher than the 
Hawley, Gizis. & Reidl (1996) average of ~ 0.5, and may in 
part be due to Malmquist bias in the flux-limited ROSAT sur- 
vey favoring the inclusion of the most X-ray luminous stars, 
as well as greater variation in the ratio because of the elapsed 
time (two decades) between the ROSAT survey and our ob- 
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FIG. 19. — X-ray luminosity normalized by the flux in the infrared K s band, 
plotted as a function of the optical-to-infrared color V — K, for stars in our 
sample which have counterparts in the ROSAT all-sky points source catalog. 
The color scheme shows the Ha equivalent width; active stars are found to 
have large X-ray flux, as expected from chromospheric activity. 
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FIG. 20. — Normalized near-UV flux as a function of the optical-to-infrared 
V — K s color. The color scheme shows the strength of the Ha equivalent 
width. Closed circles show stars identified as active based on X-ray emission, 
closed circles show stars with low or no detection in ROSAT. 



servations. 

Active stars can also be identified from ultra-vio l et ex- 
cess, as suggested in Section 3.4. Shko lnik et al.l (1201 ll 
120121) showed that GALEX UV fluxes can identify young 
M dwarfs in nearby moving groups, and can identify active 
stars to larger distances. Figure |201 shows the GALEX NUV 
to 2MASS J flux ratio, for the objects with UV detections. 
The da shed line shows the selection criteria of Sh kolnik et al.1 
(1201 ll) . Dot colors represent the EWHA index values for the 
stars, while filled circles indicate objects with logFx/Ffc > 
-2.6, i.e. stars whose X-ray flux does not identify them as 
being active. Overall, there is a good correpondence between 
the different activity diagnostics. However, there are some 
stars identified as active based on UV flux that are not iden- 
tified as such from their X-ray and/or Ha emission. Again 
this suggests that a complete identification of active M dwarfs 
in the solar vicinity may require a combination of diagnostic 
features. 

In any case our survey, which combines X-ray, UV, and 



Ha diagnostics, provides a valuable subset for identifying 
low-mass young stars in the Solar Neighborhood, and may 
potentially yield new members of young moving groups, or 
even the identification of new moving groups. The last three 
columns in Table [7] display flags for stars found to be active 
from either Ha, X-ray, or UV flux. The flag indicates activity 
by a "Y". Absence of a flag does not necessarily indicate ab- 
sence of activity: the GALEX survey does not cover the entire 
sky, and the ROSAT X-ray survey is not uniform in sensitiv- 
ity, so a non-detection in either does not necessarily indicate 
quiescence. Activity diagnostics could also be time-variable. 
Ha equivalent widt h is particular ar e know to be variable on 
various timescales dBell et ail 120 121) . In any case, there is a 
good correlation between the different diagnostics. We flag 
175 stars as active based on Ha, 42 based on X-ray emission, 
and 172 based on UV excess. Overall, 252 stars are assigned 
one or more activity flags: 19 stars have all three flags on, 99 
stars get two flags, and 137 get only one. 

8. DISTANCES AND KINEMATICS 

8.1. Spectroscopic distances 

Astrometric parallaxes are available for 631 of the stars 
in our sample, spanning the full range of colors and spec- 
tral subtypes. We combine these data with our spectroscopic 
measurements to re-evaluate photometric and spectroscopic 
distances calibrations for M dwarfs in our census. Absolute 
visual magnitudes My are calculated and are plotted against 
both V-J color and spectral subtype in Figure |2T] M dwarfs 
with signs of activity (Ha, UV, X-ray) are plotted in green, 
other stars are plotted in black. The solid red lines are the 
best-fit second order polynomials, when both active and inac- 
tive stars are used, and after elimination of 3cr outliers. The 
equations for the fits, where spT is the spectral type (K7 is -1 
and M0=0, etc) are: 

Mj = 1.194+ 1.823(V-7)-0.079(y-7) 2 (22) 
Mj = 5.680 + 0.393(SpT) + 0.040(SpT) 2 (23) 

where SpT are the spectral types, and with the least-squares 
fit performed after exclusion of 3-sigma outliers. The 1 a dis- 
persion about these relationships are ±0.61 mag for (My, V-J), 
and ±0.52 mag for (My, SpT). The smaller scatter in the spec- 
troscopic relationship suggests that spectroscopic distances 
may be marginally more accurate than the photometric ones. 
We suspect that the larger uncertainty on the photographic V 
magnitudes may be the cause. 

The most notable feature in the diagrams is that active 
stars appear to be systematically more luminous than non- 
active star s. This corroborate s the observation made previ- 
ously by lHawlev et ail (|2002), as part of the PMSU survey. 
lHawlev et alj ([2002 ) found that active stars were more lumi- 
nous by 0.48mag in a diagram of Mk against Ti05 index, used 
as proxy for spectral subtype. For stars in our census, we find 
that among stars of spectral subtype M2 and earlier, active 
stars are on average 0.46mag more luminous at a given sub- 
type than non-active stars; in the color-magnitude diagram, 
bluer stars of colors V-J < 4.0 which are active are on av- 
erage 0.47mag more luminous than non-active stars of the 
sa me color. Bot h value s agree well with the values quoted 
by lHawlev et al.l (120021) . The systematic overluminosity of 
active stars is also responsible for some of the scatter in the 
color-magnitude and spectral-type magnitude relationships. If 
we exclude active stars, the scatter about the color-magnitude 
relationship fals marginally to ±0.58, and the scatter in the 
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FIG. 21. — Absolute visual (My) and infrared (My) magnitudes for the 631 
stars in our survey for which geometric parallax measurements exist in the 
literature. Top panels: absolute magnitudes agains t V—J color, w hich follow 
the color-magnitude relationship used in Lepine & Gaidos (201 1) to estimate 
photometric distances. Bottom panel: absolute magnitudes as a function of 
spectral subtype, based on the spectral-index classification described in this 
paper. Active stars are plotted in green, and are found to be overluminous at a 
given V—J color and given spectral subtype, compared with non-active stars. 
The offset in notable for stars of earlier M subtypes (or bluer V—J colors). 

spectral-type magnitude relationship falls to ±0.49. 

The offset in absolute magnitude between active and non- 
active stars suggests that spectroscopic and photometric dis- 
tances would be more accurate for active stars in our census 
if their estimated absolute magnitudes were made 0.46mag 
brighter that suggested by Equation [23] We therefore adopt 
the following relationships to be applied only on active stars 
of subtype M2.5 and earlier: 



[Mj] 



early— active 



= 0.734+ 1.823(y-7)-0.079(y-7) 2 (24) 



\Mj\ariy-act™ = 5.220 + 0.393(S/xO + 0.040(S/?r) 2 (25) 

Again we define as "active" any star which qualifies as such 
base on any one of our criteria (Ha, UV, X-ray). There 
are several reasons that would explain why active, early- 
type stars are more luminous at a given subtype. Activ- 
ity in a n early-type M dwar f could mean that the star is 
younger (Delfo sse et al.ll9 98): early-type M dwarfs with ages 
< lOOMyr are known to b e overluminous at a gi ven color, due 
to lower surface gravity (Shkolni k et al.ll2012h . Older stars 
cou ld remain active due to interaction with a close compan- 
ion dMorgan et al.H2.Q12h . in which case the active stars would 
also appear overluminous due to this unresolved companion. 
Late-type M dwarfs, however, can remain active for long pe- 
riods of time, and would thus not require the star to be young 
or have a close companion. 

Splitting the active and inactive stars results in lowering the 
scatter of non-active stars in the subtype-magnitude relation- 
ship (to ±0.5mag). Overall, our spectroscopic distances for 
non-active M dwarfs provide a lcr uncertainty of ±26% on 
the distance. For active stars, we find a scatter of ±0.6mag, 
which suggests distance uncertainties of ±32%. Our photo- 
metric distances estimated from (V,V-J) have similar though 
perhaps slightly larger uncertainties. Photometric and spec- 
troscopic parallaxes, estimated from the above relationships 
for active and non-active stars, are listed in Table [7] 

We estimated the effect of two relevant sampling biases on 
the calibration between My and V—J color. In Eddington bias, 
photometric errors scatter more numerous, bluer, and intrinsi- 
cally brighter stars to redder apparent colors than redde r stars 
are scattered to bluer apparent colors (Eddington 191 3b . The 
net effect is to make stars at a given apparent color appear 
more luminous than they are. In Lutz-Kelker (LK) bias, a 
form of Malmquist bias, errors in trigonometric prallax will 
scatter more numerous and more distant stars with lower par- 
allax to higher apparent parall ax values, mak ing them appear 
less luminous than they are dLutz & Kelkeri [19731) . By tak- 
ing the derivative of My with respect V—J color, multiplying 
by the derivative of the number of stars in our ./-magnitude- 
limited catalog with respect to My, assuming that the errors 
in V—J are gaussian-distributed with standard deviation ay-j, 
and integrating over the distribution, we find the Eddington 
bias in My to be; 



A £ = -In 10 [1.918 = 0. 178 (V-/)] 2 (0.6 



-0.4 7 )4_ y , (26) 



where 7 is the power-law index of a luminosity function for M 
stars which we take to be 0.325. Performing a similar deriva- 
tion for the effect of L-K bias on My, we find: 



15 

inTo c 



(27) 



where av is the fractional error in parallax. Using published 
parallax errors for our Hipparcos stars and adopting a con- 
servative Oy-j = 0.05, we find that L-K bias usually domi- 
nates over Eddington bias and that 74% of our stars have a 
total bias of less than +0.2 magnitudes. A running median 
(N = 100) vs. V-J color is highest (~ 0.15) for the bluest 
(V -J = 2.7) stars and falling to less than +0.05 magnitudes 
for V—J > 3.3. An analogous analysis can be performed for 
the bias in My vs. spectral type, with a similar result. To de- 
bias values of My, these values should be subtracted from our 
calibration but we do not perform that operation here because 
of the small magnitude of the effect, which would overesti- 
mate distances by about 2% on average. The correction would 
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FIG. 22. — Top: distribution of spectroscopic distances for the stars in our 
survey, shown for three ranges of spectral subtypes. Early-type stars are 
clearly sampled over a larger volume, which explains why they dominate 
our survey, bottom: distribution of Galactic scales heights for the same stars, 
assuming that the Sun is hovering 15 pc sbove the Galactic midplane. As 
expected from our magnitude-limited sample, stars of later spectral subtypes 
(and lower luminosity) are found at shorter distances. Our survey samples a 
region well within the Galactic thin disk. 

also seem negligible compared with the intrinsic scatter in our 
color-magnitude and subtype-magnitude relationships are of 
order ±0.5 - 0.6, much larger than the L-K correction. 

Figure [22] shows the distribution of photometric dis- 
tances for our complete sample using Equation [23] and the 
My=My(V-J) color-magnitude relationship. The spectral sub- 
types are plotted in separate colors and demonstrates that the 
earlier-type stars, which are intrinsically brighter, are sampled 
to significantly larger distances compared with the later-type 
stars. In the 20pc volume, the M3-M4 stars still appear to 
dominate. We also plot the Galactic height of the stars in 
our sample, adopting a Galacti c height of 15 pc for the Sun 
(Cohenl [19951: iNg et al.1 119971; iBinnev. Gerhard. & Spergel 
1997). It is clear that our survey is largely contained within 
the Galactic thin disk, and barely extends south of the mid- 
plane. This is consistent with the relative absence of metal- 
poor stars associated with the thick disk and halo. 

8.2. Kinematic analysis 

Accurate radial velocities are not available for most of the 
stars in our sample, which prevents us from calculating the 
full (U,V,W) components of motion for each individual star. 
However, it is possible to use the distance measurements (for 
stars with parallaxes) or estimates (for stars with no paral- 
lax), and combine them to the proper motions to evaluate with 
some accuracy at least two of these components for each star. 
More specifically, we calculate the (U,V,W) components by 
assuming that the radial velocities Ry = 0. We then consider 
the (X,Y,Z) positions of the stars in the Galactic reference 
frame, from the distances and sky coordinates. For stars with 
the largest component of position in +X or -X, the radial ve- 
locity mostly contibute to U, and has minimal influence on the 
values of V and W. Likewise stars with the largest component 



of position in +Y or -Y (+Z or -Z) make good tracers of the 
velocity distribution in U and W (U and V). We use this to as- 
sign any one of (U,V) or (U,W) or (V,W) velocity component 
doublet to every M dwarf in our catalog. Estimated values 
of the components of velocity are listed in Table [7] For each 
star, one of the components is missing, which is the compo- 
nent that would most depend on the radial velocity compo- 
nent based on the coordinates of the star. Again, the other two 
components are estimated only from proper motion and dis- 
tance. For the distance, we use the trigonometric parallaxes 
whenever available; otherwise the spectroscopic distances as 
are used, based on the raltionships described in the previous 
section. 

The resulting velocity distributions are displayed in Fig- 
ure [23] We measure mean values of the velocity components 
using all allowable values and find: 

<{/ >=-8.1kms" 1 ,crf/ = 32.8kms" 1 , 

< V >= -17.0 kms" 1 , oy = 22.8 kms" 1 , 

< W >= -6.9 kms" 1 , a w = 19.3 kms" 1 . 

The values are also largely consistent with those found from 
the P MSU survey and desribed in (lHawley. Gizis. & Reidl 
1996). They are also remarkably simil ar to the moments 
of the velocity components calculated by iFuchs et ail (120091) 
for SDSS stars of the Galactic thin disk, and which are 
< U >= -8.6, cry = 32.4, < V >= -20.0,cri/ = 23.0, and < 
W >= -7 A,(Tw = 18.1. The agreement suggests that our 
distance estimates are reasonably accurate, and it corrobo- 
rates earlier results about the kinematics of the local M dwarf 
population which indicate a larger scatter of velocities in U. 
The mean values of < U > and < V > are consistent with 
the offsets from the loc al standard of rest as described in 
iDehnen & Binneyl ([1998) 

Active stars are found to have significantly smaller disper- 
sions in velocity space. All 252 stars with at least one activity 
flags (i.e. stars found to be active either from Ha, X-ray flux, 
or UV excess) are plotted in red in Figure |23] Those active 
stars have first and second moments: 

< U >= -9.3 kms" 1 , au = 25.2 kms" 1 , 

< V >= -13.4 kms" 1 , cr v = 16.8 kms" 1 , 

< w >=-6.6 kms _1 ,o-iv = 15.3 kms" 1 . 

The values are consistent with (lHawley. Gizis. & Reidll 19961) . 
who also reported that active M dwarfs tend to have a smaller 
scatter compared with inactive M dwarfs. The smaller disper- 
sion values suggest that these active stars may be significantly 
younger than the average star in the Solar Neighborhood. 

In any case, we also note that the velocioty space dis- 
tribution is non-uniform and shows evidence for substruc- 
ture. Our M dwarf data shows ve locity-space subs t ructur e 
as that observed and described in iNor dstrom e t~aT1 (12004b ; 
iHolber g, Nords trom. & Andersenl (120081) for solar-type stars 
in the vicinity of the Sun. This substructure is some- 
times referred to as "streams" or "mov ing groups", al- 
though an analysis by iBovy & Hogl (120101) shows that these 
groups do not represent coeval populations arising from 
star-formation episodes. The velocity-space substructure 
is more likely transient and associated with gravitational 
perturbatio ns which are the s ignatu re of the Galactic spi- 
ral arms Ouillen & Minchev (2005) and the Galactic bar 
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FIG. 23. — Velocity-space projections for the M dwarfs in our survey. Velocies are calculated based on photometric distances and proper motions alone (no 
radial velocities used). Each star in our census is thus displayed in only one panel, which corresponds to the projection in which the radial velocity of the star has 
the smallest contribution. Stars with significant levels of Ha emission, i.e. chromospherically active M dwarfs, are plotted in red. 

SUPERBLINK proper motion catalog used in the LG201 1 se- 
lection. This kinematic bias means that the current subset of 
M dwarfs monitored for exoplanet programs suffers from the 
same kinematic bias, which could possibly introduce age and 
metallicity selection effects. 



iMinchev et al.l (120121) . A simple description of the velocity- 
space distribution in terms of means values and dispersions, or 
as velocity ellipsoid, is therefore only a crude approximation 
of a more complex and structured distribution. 

Finally, we note that the stars of our catalog that 
were previously part of the CNS3 and stars with mea- 
sured trigonometric parallaxes (e.g. from the Hippar- 
cos catalog) tend to have larger velocity dispersions, with 
(cjf/,ov,<7w')=(38.4,25.9,23.1), while the newer stars have 
(cjf/,ov,<7w')=(26.2,19.4,15.3). The difference could be due to 
systematic underestimation of the photometric/spectroscopic 
distances, but a more likely explanation is that the CNS3 
and parallax subsample suffers from proper motion selection. 
This is because most of the CNS3 stars and M dwarfs moni- 
tored with Hipparcos were selected from historic catalogs of 
high proper motion stars, which have a higher limit than the 



9. CONCLUSIONS 

We have now compiled spectroscopic data for a nearly com- 
plete list of M dwarfs in the northern sky with apparent magni- 
tudes J < 9. Our survey identifies a total of 1,403 very bright 
M dwarfs. Our new catalog provides spectral subtypes and 
activity measurements (Ha emission) for all stars, as well as 
a rough indicator of metallicity in the guise of the £ parameter, 
which measures the ratio of TiO to CaH bandstrengths. Only 
one of the stars in the survey is unambiguously identified as a 
metal-poor M subdwarf (PM 120050+5426 = V15 13 Cyg). 

Our target stars were identified from the all-sky catalog of 
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bright M dwarfs presented in tLepine & Gaidosl (1201 lb . As 
such, our spectroscopic survey suffers from the same selection 
effects and completeness issues. The completeness and bias 
of the SUPERB LINK proper motion survey, from which these 
stars w ere selected is discussed a length in Lepi ne & Gaidosl 
(1201 II) . In the northern hemisphere, the SUPERBLINK cata- 
log is complete for proper motions fi > 40 mas" 1 . We show 
in §4 that there is a kinematic bias in the catalog which ex- 
cludes stars with very low transverse motions (in the plane 
of the sky), but the low proper motion limit means that less 
than 5% of stars within 65 parsecs of the Sun are in fact ex- 
cluded in the selection. In addition, we estimate that «5% 
of the nearby, bright M dwarfs may have escaped our target 
selection scheme due to faulty magnitudes. Therefore, we es- 
timate that our census most likely include > 90% of all ex- 
isting M dwarfs in the northern sky with J < 9. Early-type 
K7-M1 dwarfs have absolute magnitudes Mj pa 5.5, and our 
J < 9 sample thus identifies them well to a distance of about 
50pc, as confirmed in Figure [22] On the other hand, later type 
M3-M4 dwarfs have Mj w 8 and thus only those at very close 
distance range (<15 pc) will be included in the catalog. Their 
completeness will however be very high because the proper 
motion bias excludes less than 1% of the stars within that dis- 
tance range. In any case, the different survey volumes for 
early-type and late-type stars means that our survey favors the 
former over the latter by a factor of about 35 to 1 . It is thus no 
surprise that our spectroscopic catalog is dominated by early- 
type M dwarfs. 

An important result of our spectroscopic analysis is the 
identification of systematic errors in the spectral indices, 
which measure the strenghts of the CaH, TiO, and VO molec- 
ular bands. Systematic offsets between data obtained at MDM 
Observatory and at the University of Hawaii 2.2-meter tele- 
scopes, as well as offsets between these and the values mea- 
sured for the sames stars in the Palomar-MSU survey of 
iReid. Hawley. & Gizisl (Q~995), indicate that these spectral in- 
dices are susceptible to spectral resolution and spectrophoto- 
metric calibration, such that using the raw measurements may 
result in systematic errors in evaluating spectral subtypes and 



the metallicity £ parameter. In Section 3.2 we outline a proce- 
dure for calculating corrected indices, based on a calibration 
of systematic offsets between two observatories. A proper 
calibration requires that large numbers of stars be re-observed 
every time a new observatory and/or instrumental setup is 
used, in order to calibrate the offsets and correct the spec- 
tral indices. Only the corrected spectral indices can be used 
reliably in the spectral subtype and £ relationships, which are 
calibrated with respect to the corrected values. We adopt the 
Palomar-MSU measurement as our standard of reference for 
the spectral indices, and correct our MDM and UH values ac- 
cordingly. 

In the end, this catalog provides a useful list of targets for 
exoplanet searches, especially those based on the radial ve- 
locity variation method. Current methods and instruments re- 
quire relatively bright stars to be efficient, and the stars pre- 
sented in our spectroscopic catalog all constitute targets of 
choice, having been vetted for background source contami- 
nation. Our accurate spectral types will be useful to guide 
radial velocity surveys and selct stars of comparatively lower 
masses. 

We also provide diagnostics for chromospheric activity 
from Ha emission, X-ray flux excess, and UV excess. Be- 
sides being useful to identify more challenging sources for ra- 
dial velocity surveys, they also isolate the younger stars in the 
census. Follow-up radial velocity observations could tie some 
of the stars to nearby moving groups, and these objects would 
be prime targets for exoplanet searches with direct imaging 
methods. 
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